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ABSTRACT 

Using the Yale stellar evolution code, we have calculated theoretical models for nearby stars 
with planetary-mass companions in short-period nearly circular orbits: 51 Pegasi, r Bootis, v 
Andromedae, p 1 Cancri, and p Coronae Borealis. We present tables listing key stellar parameters 
such as mass, radius, age, and size of the convective envelope as a function of the observable 
parameters (luminosity, effective temperature, and metallicity) , as well as the unknown helium 
fraction. For each star we construct best models based on recently published spectroscopic 
data and the present understanding of galactic chemical evolution. We discuss our results in 
the context of planet formation theory, and, in particular, tidal dissipation effects and stellar 
metallicity enhancements. 

Subject headings: Planets and Satellites: General — Stars: Planetary Systems; Stars: individual 
— 51 Pegasi, r Bootis, v Andromedae, p 1 Cancri, p Coronae Borealis 

1. Introduction and Motivation 

The detection of planets outside the Solar System constitutes one of the most exciting recent 
developments in astronomy and astrophysics. These discoveries of extrasolar planets will lead to significant 
improvements in our understanding of many processes related to planet and star formation, as well as 
deeper questions such as the existence of extraterrestrial life in the Universe. We now know of more planets 
outside the Solar System than inside (including planets around pulsars; Wolszczan 1994). Several groups 
have reported detections of Jupiter- type planets around nearby solar-like stars (see Table 1), and new 
announcements continue to arrive every few months. We expect that radial velocity surveys already in 
progress (Marcy et al. 1997; Korzennik et al. 1997; Mayor & Queloz 1997; Cochran et al. 1997a; Butler 
et al. 1998) will discover many more in the near future. Other techniques such as photometry (Borucki 
& Summers 1984) and astrometry (Gatewood 1996; Pravdo & Shaklan 1996), including interferometric 
astrometry (Colavita & Shao 1994) and, in the future, space-based interferometric astrometry (Unwin et 
al. 1996; Boden et al. 1996), may lead to additional detections. In addition, observations of protostellar 
disks (Mannings 1998) and T Tauri stars (Beckwith & Sargent 1996) can provide important information 
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about planet formation. In the long term, space based interferometry appears very promising (Unwin et 
al. 1996; Boden et al. 1996). 

This paper is the first of a series in which we use the Yale Rotating Evolution Code (the YREC code 
developed by Demarque and collaborators; see Guenther et al. 1992) to calculate theoretical models for 
the present structure and past evolution of stars known to harbor planets, using recently measured stellar 
parameters and Hipparcos parallaxes as constraints. Our goal is to obtain a set of self-consistent theoretical 
models that will provide the best possible values (with realistic error estimates) for all derived stellar 
parameters, such as mass and radius, which are crucial for a number of theoretical investigations of these 
systems. In this first paper we present a general description of the stellar evolution code and our method 
for constructing the theoretical models. 

1.1. Tidal Dissipation 

Our models will provide a detailed description of the stellar convective zone, which is essential when 
calculating tidal dissipation through eddy viscosity (e.g., Zahn 1977; Zahn & Bouchet 1989). Tidal 
dissipation may have played a crucial role in the formation and evolution of the 51 Peg- type systems. 
Despite observational biases, the clustering of orbital periods around 3 — 4 d for these systems is probably 
not a coincidence. The San Francisco State University (SFSU) extrasolar planet search (Marcy et al. 1997) 
is capable of detecting periodic radial velocity variations with an amplitude of 12 m s _1 , corresponding to 
a lMj planet orbiting at 1 AU from a 1 Mq star given a favorable inclination. Additional factors, such 
as the variability of stellar photospheres and the longer time span of observations necessary to confirm 
the Keplerian nature of the radial velocity variations, can also inhibit detection of planets with larger 
semi-major axes. Nevertheless, the present low yield of extrasolar planets with larger semi-major axes is 
likely indicative of a real pattern. Periods of 3 — 4 d correspond to the onset of strong tidal dissipation in 
such systems. When the star is rapidly rotating, tidal torques may prevent the planet from spiraling all the 
way into the star and being destroyed (Lin, Bodenheimer, & Richardson 1996; Trilling et al. 1998). This 
may be essential in formation scenarios invoking slow inward migration of giant planets in a protostellar 
disk. Alternatively, tidal dissipation effects may circularize a highly eccentric orbit produced by a dynamical 
interaction (Rasio & Ford 1996). In addition, for a slowly rotating star (most of the stars in Table 1 are 
known to be slowly rotating today), tidal dissipation leads to (possibly rapid) orbital decay, and therefore 
the survival of the system to the present provides an important constraint on theoretical models (Rasio et 
al. 1996). In some cases this may also provide an upper limit on the companion mass (see §4.1). 

Tidal dissipation is particularly interesting yet difficult to study theoretically in these systems since the 
tidal pumping period (half the orbital period) is short compared to the typical eddy turnover timescale in 
the stellar convective envelope. In this regime, the efficiency of tidal dissipation can be drastically reduced 
(Goldreich & Keeley 1977), although the details of the theory in this case are rather controversial (see 
Goodman & Oh 1997 and references therein). The 51 Peg-type systems therefore provide important new 
constraints on the theory of tidal dissipation in a regime where it is particularly poorly understood. A 
better understanding of tidal effects (such as binary circularization and spin synchronization) in this regime 
would be applicable to a number of other important systems, including prc-main-scquence (hereafter PMS) 
binaries (Zahn & Bouchet 1989). The 51 Peg system has a very short-period (4.2 d) orbit, which makes it 
particularly interesting in the context of tidal dissipation theory (Rasio et al. 1996; Marcy et al. 1997). 
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1.2. Planet Formation 

The properties of the new planets are surprising (see Table 1). Most are Jupiter-mass objects in very 
tight circular orbits or in wider eccentric orbits. The standard model for planet formation in our Solar 
System cannot explain their orbital properties (see, e.g., Lissauer 1993; Boss 1995). According to this 
standard model, planetary orbits should be nearly circular, and giant planets should be found at large 
distances (^ 1 AU) from the central star, where the temperature in the protostellar nebula is low enough 
for icy materials to condense (Boss 1995, 1996). These simple predictions of the standard model for the 
formation of the Solar System are at odds with the observed parameters of most detected extrasolar planets. 
With only one exception (47 UMa), the new planets all come within 1 AU of the central star. Three 
planets (51 Peg, r Boo, v And) are in extremely tight circular orbits with periods of only a few days. Two 
planets (p 1 Cnc and p CrB) have circular orbits with somewhat longer periods, of order tens of days. Three 
companions with wider orbits (16 Cyg B, 70 Vir, and HD 114762) have very large eccentricities (~ 0.5). A 
number of different theoretical scenarios have been proposed to explain the unexpected orbital properties 
of these extrasolar planets. Our calculations will provide accurate stellar parameters that can be used to 
obtain better constraints on these formation scenarios. 

Four specific mechanisms have been suggested as possible processes to bring a giant planet into a 
short-period orbit around its star. One mechanism is a secular interaction with a distant binary companion. 
If the orbit of a wide binary is inclined relative to a planet's orbit by more than ~ 40°, the relative 
inclination of the binary star can couple to a secular increase in the eccentricity in the planet's orbit 
(Holman, Touma, & Tremaine 1996; Mazeh, Krymolowski, & Rosenfeld 1996). The amplitude of the 
eccentricity perturbation depends on the relative inclination of the orbits, but is independent of the mass 
of the binary companion. If the orbital planes have a very high relative inclination, then an eccentricity 
approaching unity can be induced. In some cases this may cause a collision with the star. However, if 
dissipation is significant (in the star, in the planet, or in a disk), then the orbit could circularize at a small 
distance. Similarly, if the primary star has a significant quadrupole moment, then tidal dissipation in the 
star could stop the growth of the eccentricity oscillations and drive a gradual decrease in the semi-major 
axis of the planet's orbit (Kiseleva & Eggleton 1997; Eggleton & Kiseleva 1997). This mechanism of secular 
perturbations from a distant binary companion can also produce significant eccentricities for the longer 
period extrasolar planets. It is interesting to note that four of the five planets with semi-major axes, 
a <, 0.25 AU arc in wide binary systems (Hoffleit & Warren 1991). 51 Peg has been searched extensively 
for a binary companion (Marcy et al. 1997), but none has been found. If the relative inclination less than 
~ 40°, then the amplitude of the eccentricity oscillations becomes sensitive to the mass of the perturbing 
body, the ratio of the semi-major axes, and the eccentricity of the outer orbit. For a binary companion 
of stellar mass, only a small range of parameters will produce large eccentricities without disrupting the 
system. The large mass ratio of either star to the planet makes it easy for the system to be disrupted. 
Since the planet is the least massive of the three bodies by many orders of magnitude, the planet is almost 
certainly the body to be ejected, if the system dissociates. 

Dissipation in the protostellar nebula is a second possible mechanism for forming a short period planet. 
Since the orbital migration of the planet would tend to accelerate with decreasing separation from the star, 
the dissipation has to switch off at a critical moment for a planet to end up so close to its parent star without 
being disrupted. Possible mechanisms for stopping the inward migration include Roche lobe overflow and 
tidal coupling to a rapidly rotating star (Trilling et al. 1998; Lin et al. 1996). Another possibility is that 
the migration stops when the planet arrives at the inner edge of a disk limited by a magnetosphere around 
the star (Lin et al. 1996). 
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A resonant interaction with a disk of planetesimals is another possible source of orbital migration, but 
this requires a very large protoplanetary disk mass if a ~ lMj planet is to migrate inwards all the way to 
~ 0.1 AU (Murray et al. 1998). The advantage, however, is that the migration is halted naturally at short 
distances when the majority of perturbed planetesimals collide with the star rather than escaping on nearly 
parabolic orbits. Wide eccentric orbits can also be produced for planets more massive than ~ 3Mj. It 
seems that such a massive disk would be likely to produce more than one planet, and thus it is important 
to understand how a second or third planet would affect this scenario. 

The fourth mechanism is based on dynamical instabilities in a system originally containing multiple 
giant planets of comparable masses (Rasio & Ford 1996). If either the orbital radii evolve secularly at 
different rates (significant orbital migration is thought to have occurred in the outer Solar System; see 
Goldreich & Tremaine 1980; Malhotra 1995) or if the masses increase significantly as the planets accrete 
their gaseous envelopes (Lissauer 1993), then the orbits could become unstable. Alternatively, if either 
the second or third mechanism for producing short period planets were truncated when the disk mass 
dwindled, they could leave two planets in dynamically unstable orbits. Similarly, a wide binary companion 
could drive a secular increase in the planets' eccentricity until they became dynamically unstable. In any 
case, the evolution can lead to a dynamical instability of the orbits and a strong gravitational interaction 
between two planets (Gladman 1993; Chambers, Wctherill, & Boss 1996). This interaction can lead to the 
ejection of one planet, leaving the other in an eccentric orbit. If the pericenter distance of the inner planet 
is sufficiently small, its orbit can later circularize at an orbital separation of a few stellar radii (Rasio et al. 
1996). This mechanism can produce eccentric systems in two different ways. After many encounters, one 
planet can be ejected from the system leaving the other planet in an eccentric orbit by itself (Rasio & Ford 
1996; Katz 1997). In this scenario it should be expected that the second planet has been ejected and is 
no longer in the system (Black 1997). Alternatively, if the two planets collide, they can produce a more 
massive planet in an eccentric orbit (Rasio & Ford 1996; Lin & Ida 1997). Dynamical instabilities in systems 
containing more than two giant planets of comparable masses have also been studied. Wcidcnschilling & 
Marzari (1996) have obtained numerical results for the case of 3 planets, while Lin & Ida (1997) performed 
simulations for systems containing up to 9 planets. In this case successive mergers between two or more 
planets can lead to the formation of a fairly massive (<; 10 Mj) object in a wide, eccentric orbit. While it is 
almost certain that this mechanism operates in many systems with multiple planets, it is not clear that it 
can produce a fraction of 51 Peg-like systems as large as observed. Extensive numerical simulations will be 
necessary to obtain good statistics on this theory (Ford & Rasio 1998). Observational selection effects must 
also be better understood for a meaningful comparison with the properties of the detected systems. 

All these dynamical processes can also affect the evolution of the central star. Marcy et al. (1997) and 
Drake et al. (1998) argue that the observed rapid rotation of r Boo is likely caused by the tidal interaction 
between the star and its close planetary companion (our results do not support this interpretation; see §4.1). 
Several stars with short-period planets have high metallicities ([7^] lb 0.2; see Gonzalez 1997a, 1998ab 
and Table 2). Gonzalez (1998a) proposes that their metallicities have been enhanced by the accretion of 
high-Z material. Alternatively, the correlation could arise, if metal-rich stars have metal-rich disks which 
are more likely to form planets. Thus, understanding the relationship between metallicity and the existence 
of short-period planets may be important in constraining the mechanisms which produce these planets. We 
discuss this in more detail after constructing models for these stars (See §4.2). 



2. Constructing Stellar Models 
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2.1. The Code 

We use the Yale Rotating Evolution Code in its non-rotating mode to calculate stellar models. YREC 
is a Henyey code which solves the equations of stellar structure in one dimension. The chemical composition 
of each shell is updated separately using the nuclear reaction rates of Bahcall & Pinsonneault (1992). The 
initial chemical mixture is the solar mixture of Anders & Grevesse (1989), scaled to match the metallicity 
of the star being modeled. For regions of the star where logT(K) > 6, all atoms are assumed to be fully 
ionized. For regions where logT(K) < 5.5, particle densities are determined by solving the Saha equation 
for the single ionization state of hydrogen and the metals, and for the single and double ionization states 
of helium. In the transition region between these two temperatures, both formulations are weighted with a 
ramp function and averaged. The equation of state includes both radiation pressure and electron degeneracy 
pressure. We use the latest OPAL opacities (Iglasias & Rogers 1996) for the interior of the star down to 
temperatures of logT(K) — 4. For lower temperatures, we use the low-temperature opacities of Alexander 
& Ferguson (1994). We use an Eddington T-r relationship for the atmosphere, and, where appropriate, the 
standard Bohm-Vitense mixing length theory of convection (Cox 1968; Bohm-Vitense 1958) with the ratio 
of the mixing length to the pressure scale height, a=1.70. This value of a, as well as the solar hydrogen 
abundance, Xq — 0.697, was obtained by calibrating models against observations of the present day Sun 
using the aforementioned physics. 

Astroseismology has made it possible to test models of stellar interiors directly. Models of the Sun 
(see, e.g., Guenthcr ct al. 1992) are constantly being refined with newer and more sophisticated input 
physics, including diffusion of helium and heavy elements (Guenther & Demarque 1997) and detailed 
hydrodynamic calculations of the surface convection zone (Demarque et al. 1997). YREC has also been 
used to calculate models of other single field stars. Most recently, observations of p-mode oscillations in 77 
Bootis, a nearby subgiant, prompted detailed stellar modeling of that star (Guenther & Demarque 1996). 
This work forms the basis of our method. The modeled frequency spectra of the Sun and rj Bootis agree 
quite well with observations. YREC has also been used to model Procyon (Guenther & Demarque 1993) 
and both components of the binary star a Centauri (Edmonds et al. 1992). Although stellar models have 
also predicted seismology spectra for Procyon and a Centauri, no convincing observations of the frequency 
spectra have yet been made. 

In constructing stellar models a variety of formalisms are available for the treatment of the equation of 
state, the stellar atmosphere, and convection. We adopted formalisms which result in a good solar model, 
and are standard enough to be applicable to all low-mass stars. We have neglected non-standard effects 
such as diffusion of helium and heavy elements, rotation, and magnetic fields since these processes do not 
significantly change the overall properties of low-mass stars like the sun. Solar models constructed by 
YREC with the set of standard formalisms adopted here reproduce the observed solar p-mode spectrum to 
within 1% (Guenther ct al. 1992). The incorporation of newer physics can reduce this difference by a few 
tenths of a percent. 

We have chosen to use the Saha equation of state rather than MHD (Mihalas et al. 1988) or OPAL 
(Rogers & Iglesias 1994). The MHD and OPAL equations of state are very similar to each other, and differ 
from the YREC implementation of the Saha equation by less than 1%. The OPAL equation of state results 
in a slightly deeper convection zone for the sun, and a worse fit to the observed p-mode spectrum (Guenther 
et al. 1996) compared to the Saha equation of state. 

For the treatment of convective zones we have adopted the standard Bohm-Vitense mixing length 
theory. Some other theories, such as that of Canuto & Mazzitelli (1992), produce solar models which have 
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p-mode spectra in slightly better agreement with observations to the observed spectrum. Currently, a 
number of groups are using 3-D numerical simulations of turbulent convection to develop a more realistic 
description of stellar convective zones (e.g. Kim & Chan 1998). However, the results of these simulations 
are not yet available in sufficient detail to be incorporated into general stellar models. 

The mixing length is chosen so that our solar model reproduces the solar luminosity and solar radius 
at the solar age. Our best value of the mixing length parameter (ratio of mixing length to the local pressure 
scale height) is 1.7. To test the effect of the mixing length on the stellar parameters, for several of our best 
models we varied the mixing length from 1.4 to 2.0 (a large range). We found that the luminosity of the star 
typically varied by ~ 3%, the radius by ~ 5%, the mass of the convective envelope by ~ 35%, and the radius 
of the convective envelope by ~ 2%. We have also performed a number of tests to check that the choice of 
atmosphere model does not affect our results significantly. For examples, we find that using the Kurucz 
tabulated atmospheres (Kurucz 1991) rather than the Eddington T-r relation increases the calculated 
effective temperatures and increases the calculated radii by less than 1%. We can test our complete set of 
physical parameters and assumptions by comparing the calculated depth of the solar convection zone to 
that derived from the observed solar p-mode spectrum. Christensen-Dalsgaard et al. (1991) give the radius 
at the base of the solar convection zone to be 0.713 ±0.003 R Q . Our calibrated model produces a convection 
zone which begins at 0.647 R Q , a difference of ~ 10%. In summary, we expect that the calculated stellar 
radii and temperatures will be accurate to within about 5%, and the calculated size of the convection zones 
will be accurate to within 20%, including all possible changes to the physics of the models. 



2.2. Iteration Method 

We constructed models of O.6M , 0.8M©, l.OM , 1.2M , 1.4M , and 1.6M PMS stars by solving 
the Lane-Emden equation for a polytrope of index n = 1.5. YREC evolved the models to the zero-age 
main sequence (hereafter ZAMS). These ZAMS models serve as the starting point for all our subsequent 
modeling. 

For each model we construct, we first specify a mass and metallicity and then evolve a ZAMS model 
in an attempt to match a given effective temperature T e g and luminosity L. First, YREC chooses a ZAMS 
model with a similar mass and metallicity and scales the ZAMS model to match the desired values of M* 
and Z (indirectly affecting X and Y). The scaled ZAMS model is numerically relaxed before YREC begins 
to evolve the model to the desired radius (determined from effective temperature and luminosity). YREC 
iterates the above procedure, returning to the ZAMS model, but rescaling to new values X and Y so as to 
improve the match to the desired luminosity and radius. 

Thus, a given run fits a model to the desired luminosity and radius (and hence temperature). For such 
a run we hold M* and Z constant, explicitly vary Y, and let X, the age, and the other stellar parameters 
vary as a result of the changing composition. This procedure is repeated for several masses yielding a set 
of models which share the same specified metallicity, luminosity, and temperature, but differ in mass, age, 
composition, and other parameters. 

Since we vary the helium content in our models, it is possible to obtain models that match the observed 
parameters but have unrealistic compositions. To select our final "best" ' model, we consider the ratio 



SY _ Y-Yq 



5Z Z-Zq 

where Y and Z are the helium and heavy element abundances of our model. Based on studies of galactic 
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evolution, including H II regions (Pagcl et al. 1992) and low-metallicity blue compact galaxies (Izotov et 
al. 1997), we impose the canonical constraint 4^ = 2.5 ± 1 (Bressan et al. 1994; Edvardsson et al. 1993). 
For our "best" model, we interpolate to find a mass and age corresponding to |^ = 2.5 for the specified 
parameters. Similarly, we present models for |^ equal to 1.5 and 3.5 from which we determine theoretical 
error bars given a set of assumed observational parameters. Due to the physical scatter and the possibility 
of systematic errors in the determination of |^-, we also include models for |^ equal to 5.0 and 0.0. The 
latter may also be useful for comparisons with studies in which Y — F Q is assumed. Fortunately, the 
theoretical uncertainties derived from the variation of |^ are, in general, smaller than the uncertainties 
in the stellar observations, as will be discussed in the next section. The above procedure is repeated for 
each set of observed parameters, Z, T e g, and L. For some sets of observed parameters, YREC is unable to 
converge on a single self-consistent model. 



2.3. Observational Data 

The above procedure for constructing a set of models to match a real star requires three basic input 
parameters from observations. The luminosity, L, effective temperature, T c g, and metallicity, Z, must be 
known accurately. The observational uncertainty in these quantities limits the accuracy with which we can 
derive the values of other stellar parameters. As more observational data become available, these can be 
used to constrain our models more tightly. 



2.3.1. Luminosities 

To minimize systematic errors, we consistently use the stellar luminosities and parallaxes obtained 
from the Hipparcos catalog. The recent release of the Hipparcos data (ESA 1997, Perryman et al. 1997) is 
extremely useful for this type of study, since accurate distances to all the relevant stars are now available. 
The Hipparcos parallax data have greatly improved astronomical distances and tightly constrained 
luminosities, nearly removing an entire degree of freedom from the models. Since only the brightest (and 
therefore the closest) solar-type stars are targets for present spectroscopic surveys, all the solar-type stars 
presently known to have planets are nearby and have distances directly determined by Hipparcos parallax 
data. This distance is combined with the Hipparcos apparent visual magnitude to obtain a visual luminosity. 

One minor complication is the application of bolomctric corrections to convert visual luminosities 
to integrated luminosities. We estimate bolometric corrections by interpolating in [Fe/H], logg, and T c ff 
across color calibration grids used for the Yale Isochrones (Green et al. 1987). This method of estimating 
bolomctric corrections does introduce a dependence of our results on log g. However, bolometric corrections 
depend only weakly on log<?, and the uncertainty in logg hardly introduces any uncertainty in the 
bolometric correction. Comparing our bolometric corrections with other tables indicates that there is a 
~ 5% systematic uncertainty in our luminosities due to the potential systematic error in the bolometric 
corrections. From our full set of models we see that the uncertainty in the integrated luminosity is normally 
insignificant compared to other sources of error. 
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2.3.2. Spectroscopic Data 

In addition to the luminosity, our models also requires a knowledge of each star's metallicity and 
effective temperature. Effective temperatures and metallicities can be obtained from either photometric 
or spectroscopic data, although the latter are normally more accurate. In this paper we have chosen 
values for the necessary parameters from the literature. To minimize systematic errors, we decided to use 
temperatures, metallicities, and surface gravities obtained from a single source, the recent spectroscopic 
observations of all the solar-like stars with planets by Gonzalez (1997ab, 1998ab). His determinations of 
T ff, loggf, [Fe/H], and the microturbulence parameter ( t , are from self-consistent iterative solutions which 
matched high resolution spectra of Fe I and Fe II lines to Kurucz (1993) model atmospheres. Gonzalez 
(1998) quotes typical errors of 75 K, 0.06, and 0.05 for T c s, [Fe/H], and iogg, respectively. These formal 
errors may be somewhat smaller than the actual uncertainty in the values, but the general agreement with 
other observations is reassuring. We construct additional models when other observations significantly 
deviate from the Gonzalez results, or when we cannot construct models simultaneously matching all the 
observed parameters. We summarize the observed parameters of stars with planets in Table 2. We will 
list other determinations of the stellar parameters, T e s, [Fe/H], \ogg, and vsini, as we discuss each star 
individually. 

2.3.3. Other Observational Parameters 

Several other stellar parameters can be measured observationally. For example, logg is routinely 
measured spectroscopically. While determinations of logg serve as a useful tool for purposes such as 
determining bolometric corrections (see §2.3.1), present determinations (ci ogg ( cgs ) ~ 0.1) provide only loose 
constraints (~ 25%) on M*/R%. Thus, we use measurements of logc/ merely as a consistency check for our 
models. 

Analysis of Ca II H and K emission can be used to detect the stellar rotation period. When a rotation 
period is not detected, it can be predicted using an empirical relation between the Ca II flux and the 
rotation period (Noyes et al. 1984). The rotation period can be combined with the radius to yield an 
equatorial velocity, v cq . This can be compared with the observed vsini to yield the inclination angle 
between the star's equator and the line of sight. If the angle between the star's equator and the planet's 
orbital plane is presumed small, one can determine the planet's actual mass, to, from its minimum mass, 
to sin i . Unfortunately, vsini is very difficult to determine observationally, as it is delicately coupled to 
the macroscopic turbulence parameter. We list only recent determinations of vsini and take even these 
somewhat cautiously. Empirically, both the level of Ca II H and K activity and the stellar rotation 
frequency have been found to decrease with age and this relationship can be used to estimate the age of 
the star (Baliunas et al. 1995; Sodcrblom ct al. 1991). We will compare our determinations of the age with 
those predicted by the activity-age relation. 

Photometric determinations of the angular diameter are possible with the Barnes and Evans relationship 
(Barnes, Evans, & Moffet 1978; Moffett & Barnes 1979) or the more recent infrared flux method (Blackwell 
et al. 1990). Combined with parallax measurements, these yield a stellar radius. Given the relatively high 
uncertainties associated with these methods, we also use these observations only as a consistency check. 
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3. Results 

For each star we have computed a grid of models surrounding the observed values of L, T c g, and [Fe/H]. 
As more observations become available, it is hoped that our large grids of models will allow observers to 
translate their values of the observed parameters into physical parameters such as mass, radius, age, and 
size of the convective zone. In general, interpolating across grids of stellar models can be inaccurate as 
non-linearities become significant. However, we have calculated a fine enough grid of models for each star 
such that interpolation should yield reasonable results. The set of models printed here is only a small subset 
of our large grid of models. Here we present models corresponding the adopted values of T c ff, [Fe/H], and L 
and models corresponding to a 1 — a uncertainty in any one of these parameters. For each model we include 
the mass M*, the age, the mass of the convective envelope M co , the radius of the convective envelope R ce , 
the pressure scale height at the base of the convective envelope PSH, the radius i?*, and the surface gravity 
logy. The full set of models will be made available electronically. 

We summarize our results in Table 3 with determinations of the mass, radius, age, size of convective 
envelope, and the eddy turnover timescale at the base of the convective envelope for each star (See §4.1). 
The error bars indicate the range of values obtained in models with |^ £ [1.5,3.5] and each of the input 
parameters, L„, T ff, and [Fe/H] varied within the ranges specified below for each star. These error bars 
do not include the uncertainty in the choice of physical models discussed in §2.1. For our best model of 
each star we also list the dimensionless gyration radius of the star, the dimcnsionless gyration radius of the 
convective envelope (See §4.1) and the mass of the convective envelope at the ZAMS. 

3.1. 51 Pegasi 

The planet around 51 Pegasi was discovered by Mayor & Queloz (1995) and has since been confirmed 
by both Marcy et al. (1997) and Horner et al. (1997). The SFSU team has the cleanest data, giving an 
rms scatter of 5.2 m s _1 about a Keplerian fit with a semi-amplitude of 55.9 ± 0.8 m s _1 . They calculate 
an orbital period of 4.2311 ± 0.0005 d and an eccentricity of 0.012 ± 0.010. The possibility of a second 
companion has been carefully examined. The SFSU planet search should already be able to detect planets 
with msini ~ lMj within 2 AU. A longer temporal baseline of observations is necessary to extend this limit 
(Marcy et al. 1997). Several teams have searched for periodic spectral line bisector variations. Although 
there is no convincing evidence (Horner et al. 1997; Hatzes, Cochran, & Cohns-Krull 1997; Gray 1998; 
Hatzes, Cochran, & Bakker 1998), the possibility of exciting non-radial stellar pulsations in 51 Peg-like 
systems remains interesting (Willems et al. 1997; Terquem et al. 1998), and these may become detectable 
in the future. 

We calculated our models with the values from Gonzalez (1998a): T c g = 5750 ± 75 K and 
[Fe/H]= 0.21 ± 0.06. The consistent determinations of these parameters from many observations (See Table 
3) is comforting. We obtain a best model of 51 Peg with a mass of l.O5±°;^M , radius of 1.16 ± O.O7i? , 
and an age of 7.6+g'5 Gyr. We find that 51 Peg has a convective envelope (M ce ~ 0.023+QQQgM Q ) only 
slightly larger than that of the sun (M ce ~ 0.0174M Q ). The eddy turnover timescale calculated as in Rasio 
et al. (1996) is r ce ~ 18.6 ± 2.5 d which is close to the expected value for main-sequence solar-type stars, 
T ce ~ 20 d. In addition, we can follow the history of the convective envelope. We find that the mass of 
convective envelope on the ZAMS was M co ,zams — 0.037M Q . 

As a consistency check, we compare our predicted value of \ogg (cgs) = 4.33 ± 0.09 with those 
determined observationally. We find general agreement, especially when we consider that several of the 
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lower observed values have already been criticized in the literature. Fuhrmann, Pfeiffer, & Bcrnkopf 
(1997) suggest that the low value of the Gratton, Carrcnton, & Castelli (1996) study is a consequence of 
their low T e ff, and that the McWilliam (1990) result should be discounted, as it depended on an earlier 
misclassification of 51 Peg as a subgiant. Indeed, the reclassification of 51 Peg as a dwarf has been verified 
by multiple observations. Fuhrmann et al. (1997) also discount the Xu (1991) value, as it is based on low 
resolution spectra. Finally, the Edvardsson et al. (1993a, b) result is superseded by their own more recent 
observations which appear in Tomkin et al. (1997). The remaining determinations are all consistent with 
our models. 

The activity- age relation predicts an age of 10 Gyr consistent with our result for 51 Peg (Baliunas et 
al. 1997), and a rotation period of 29.6 d, while a rotation period near 37 d has been observed (Baliunas 
et al. 1996). Combining a rotation period of 34 ± 4 d with our radius, we compute an equatorial velocity, 
v cq — 1.78 ± 0.23 km s _1 . This can be combined with observed values of vsini to provide estimates of 
the inclination and hence the actual companion mass. The determination of v sini = 1.4 ± 0.3 km s _1 by 
Gonzalez (1998) is quite consistent with our calculated v cq , implying sini = 0.8 ±0.2 and m ~ 0.59±0.15Mj. 
However, most recent determinations of vsini are larger than the above computed v cq (See Table 4). It 
should be noted that both the values of Hatzes et al. (1997) and Francois et al. (1996) have already been 
corrected according to Gonzalez (1998). These measurements are either inconsistent or barely consistent 
with the radius and period, suggesting that sini ~ 1 and thus the companion mass is not very different 
from the minimum companion mass. 

3.2. t Bootis 

The SFSU planet search discovered a planet around r Boo in a 3.3128 ± 0.0002 d near-circular orbit. 
There is an rms scatter of 13.9 m s _1 with occasional episodes of greater scatter about the Keplcrian fit 
with semi-amplitude 469 ± 5 m s . The scatter is significantly above the instrumental error, but cannot be 
explained by the presence of a second planet (Butler & Marcy 1996). Several teams have searched r Boo 
for periodic spectral line bisector variations, but none have been found (See §3.1). 

We started with models based on the values of T c g and [Fe/H] from the recent spectroscopic work of 
Gonzalez (1997a). In the process, we found that our models favored lower values of T e s and [Fe/H], and 
so we considered other determinations of T c g and [Fe/H]. Upon rcanalysis of his own data, Gonzalez has 
refined his estimate of T e g to 6550 ± 100 K (Gonzalez 1997b). He claims that his data might be consistent 
with T e s = 6400 K, but that they are inconsistent with T e g — 6300 K. There may be a significant difference 
between spectroscopic and photometric determinations of T c ff and [Fe/H] (Gonzalez 1997a). The small 
number of metal-rich stars used in the calibration of photometric estimators is one possible explanation 
for the difference in spectroscopic and photometric estimates for r Boo (Gonzalez 1997a). Alternatively, if 
r Boo has a convective envelope, the tidal torque from a planet in a 3.3 d orbit could drive the envelope 
to rotate much more rapidly than normal, thereby driving activity on the stellar surface and altering the 
observed spectral characteristics of the star. However there is a more likely explanation: the high X-ray 
luminosity, radial velocity noise, young age, and observed rotation period are all consistent will r Boo 
being a young star that is still rotating rapidly. A rotation period as short as 3.3 d can cause significant 
discrepancies between spectroscopic and photometric determinations of T e g and [Fe/H]. The rapid rotation 
will broaden spectral lines and slightly redden the color. This, like the high mctallicity of r Boo, can affect 
the continuum and thus all the line depths and widths. Photometric estimates for such rapidly rotating 
stars can also be adversely affected. 
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We explored a range of T e g and [Fe/H] significantly larger than the uncertainties in the Gonzalez 
(1997a, b) observations require. YREC was unable to construct any self-consistent models of r Boo in 
the range of 1.0 to 1.8 M Q using T ef j > 6600 K and [Fe/H]> 0.13, the spectroscopic value from Gonzalez 
(1997a). If we adopt T off = 6400 ± 100 and [Fe/H]= 0.25 ± 0.09, then we arrive at a mass of 1.37 ± O.O8M , 
a radius of 1.41+QQg_R Q and an age of 1.2+Qg Gyr. Our models predict logg = 4.27to'o7' which is in 
agreement with most of the observations, but slightly inconsistent with the Gonzalez (1997a) value of 
4.5 ± 0.15. The most interesting result from our models of r Boo is the likely absence of a convective 
envelope. We find thin convective envelopes only in models with T c g <^ 6350 K or [Fc/H]<^ 0.07 (See Table 
7). 

The activity-age relation predicts an age of ~ 2 Gyr (Baliunas ct al. 1997), consistent with our age. The 
Mount Wilson HK Project has revealed that r Boo has chromospheric emission periods of approximately 
3.5 ± 0.5 d, 117 d, and 11.6 yr (Baliunas ct al. 1997). The shortest period is believed to be the rotation 
period, which corresponds very nearly to the orbital period of the planet. That has led Marcy et al. (1997) 
and Drake et al. (1998) to suggest the possibility that the star may have been tidally spun up. If we assume 
the star's rotation period is synchronized with the planet's orbital period and combine the 3.3 d rotation 
period with our determinations of radius, then we calculate v cq = 21.7+\'\ km s _1 . This is consistent with 
all determinations of wsinz and suggests that sini ~ 0.67^061 implying m = 7.1 ± 0.8Mj. 

3.3. v Andromedae 

The SFSU team detected a companion to v And in a 4.611 ± 0.005 d near-circular orbit. There is a 
12.1 m s _1 rms scatter about the Keplerian fit of semi-amplitude 74.1 ± 0.4 m s _1 . This scatter is well 
above the instrumental error and could be either intrinsic to the star or indicative of a second planet. Early 
observations reported a significant eccentricity of 0.109 ± 0.040 (Butler & Marcy 1996), but the long-term 
trend in the residuals complicates the measurement of the eccentricity. Thus, the current uncertainty is 
likely larger than the quoted error bar (Marcy 1998). 

Based on the observations of Gonzalez (1997a), we calculate a mass of 1. 34^^12-^07 radius of 
1.56to'iJii0, and age of 2.6±f Gyr. The convective envelope is very shallow with M ce ~ 0.002t°;°g2^o 
and the eddy turnover time is correspondingly short, r ce ~ 6.8tg;| ! d. 

The activity-age relation predicts an age of 5 Gyr, consistent with our results (Baliunas et al. 1997). 
Based on Ca II emission, v And is expected to have a rotation period of ~ 12 d (Baliunas et al. 1997). If 
accurate, one can infer an equatorial velocity v cq ~ 6.6 km s _1 . This value is low compared to observed 
values of vsini ~ 9 km s -1 (See Table 8). This suggests a significant error in either the observed value of 
wsini or the estimated rotation period. If we were to take the reported usini at face value, then we would 
expect a rotation period <> 6 d. 

Our models give logg = 4.18t;°Q, which is consistent with all but one of the observed values of logg. 
Blackwell et al. (1990) have used the infrared flux method to calculate an angular diameter of 1.103 ± 0.044 
mas, which, combined with the Hipparcos parallax, gives i?* = 1.60 ± 0.08i? Q , in agreement with our 
models. This corrects some previous estimates based on the Blackwell et al. (1990) data and older parallax 
data (most notably 56.8 ±4.1 mas from van Altena et al. 1995 which led to incorrect estimates of T c g and 
the radius). 
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3.4. p 1 Cancri 

The SFSU team has discovered a companion around p 1 Cnc in a near-circular orbit of period 
14.648 ± 0.0009 d. They find a Keplerian fit to the radial velocity with a semi-amplitude of 77.1 ± 0.9 m 
s _1 and a rms scatter of 12.0 m s . In addition, there appears to be a long-term trend in the residuals, 
possibly indicative of a second planet with mass ~ 10Mj and a period of ~ 20 yr (Marcy & Butler 1998). 
Previously, McAlister et al. (1993) had searched for a luminous companion in a similar orbit using speckle 
observations, but found none. 

Initially, we attempted to construct models based on the stellar parameters measured by Gonzalez 
(1998ab). However, we found that this value of the temperature, T c g = 5150 ± 75 K, is too low to match 
any model of p 1 Cnc, even when evolved for 20 Gyr. 

There are several possible explanations for the apparent inconsistency of the various available data for 
p 1 Cnc. One possibility is that p 1 Cnc is actually a subgiant. Indeed, p 1 Cnc's spectra closely matches that 
of 5 Eri, a subgiant with [Fe/H]~ —0.15 (Baliunas et al. 1997)0. Gonzalez (1998a) agrees that the spectrum 
is suggestive of a subgiant, as is the low surface gravity. However our models, as well as the isochrones used 
by Gonzalez (1998a), would predict an extremely large age, of J> 12 Gyr. 

Gonzalez (1998) proposes another possible explanation: the p 1 Cnc system may be an unresolved 
stellar binary viewed nearly face-on, despite the a priori low probability of such a viewing angle. Gonzalez 
(1998) suggests monitoring the line profile variations to test this hypothesis. We discuss this possibility 
further in §4.1.2, and find this hypothesis unattractive. 

Another possible explanation is an error in the stellar models. The observed value of [Fe/H], 
+0.45 ± 0.03 for p 1 Cnc is quite extreme and the use of 5Y/SZ = 2.5 ± 1.0 to determine Y for our models (as 
well as the isochrones of Schaller et al. (1992) and Scaerer et al. (1993a, b) used by Gonzalez 1998) may not 
be appropriate for such high metallicity stars. While this constraint on SY/SZ is believed to be reasonable 
for most solar-like stars, the calibration is based on lower metallicities and this linear relationship may not 
be adequate for high values of Z. 

Finally, it is possible that the observations have larger errors than those quoted. In particular, a higher 
T e g would restore consistency. However, since the observed value of [Fe/H] is correlated with T e g, the 
observed metallicity of p 1 Cnc would then increase from its (already high) quoted value (Gonzalez 1996, 
1998ab). If we adopt T off = 5300 ± 75 and [Fe/H]= 0.45 ± 0.03, then we find a mass of 0.95t ;JgM Q and 
a radius of 0.93^_q^Rq. The convective envelope (M ce ~ 0.046^0 006-^©) is significantly larger than solar 
and the eddy turnover time is slightly longer, t c ~ 26.7+^2 d- The age is very sensitive to Y and thus is not 
well constrained. This explanation of our difficulties in modeling p 1 Cnc is supported by the more recent 
observations of p 1 Cnc by Fuhrmann et al. (1998) and Gonzalez (1998b). Given these recent observations 
and the difficulties with the alternative explanations, we find this explanation the most attractive. Although 
the Fuhrmann et al. (1998) and Gonzalez (1998b) observations consistently support a larger T e ff, the 
unknown helium fraction still prevent us from obtaining an accurate determination of the age. Fuhrmann 
et al. (1998) obtained a maximum age based on their own observations. In addition to their observations 
having the largest T G ff, they have assumed a solar helium abundance. If we use their observational data and 
assume a solar helium abundance, then we obtain a similar upper limit on the age. However, if we use their 
observational data, but assume SY/5Z = 2.5 ± 1.0, then we can no longer impose such a constraint. Thus, 



2 Cayrel de Strobel et al. 1997 lists many values of [Fe/H] for 6 Eri, ranging from -0.27 to 0.33. 
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their determination of a maximum age for p 1 Cnc is only valid if p 1 Cnc has a near solar helium abundance. 

In light of the problems modeling p 1 Cnc, we consider several other pieces of observational data. Our 
models give log g = 4.501qq7, significantly higher than the Gonzalez (1996, 1998) value of 4.15 ± 0.05, but 
in agreement with other observations (See Table 10). This is expected since log g is correlated with T e g. 
The activity-age relation suggests an age of ~ 5 Gyr (Baliunas et al. 1997). The observed rotation period 
of 41.7 d (Baliunas et al. 1997) can be combined with our radius to obtain v cq — 1.12 ± 0.03 km s _1 , barely 
consistent with the observed value of vs'mi, 1.4 ± 0.5 km s _1 . The Barnes and Evans relationship (Barnes 
et al. 1978) gives an angular diameter of 0.790 ± 0.032 mas, which we combine with the Hipparcos parallax 
to calculate a radius of 1.06 ± 0.04i? Q , larger than the radii we find in our models. 

3.5. p Coronae Borealis 

The Advanced Fiber Optic Echelle (AFOE) spectrograph team discovered a companion orbiting p CrB 
in a 39.645 ± 0.088 d near-circular orbit. They fit a Keplerian curve of semi-amplitude 67.4 ± 2.2 m s _1 to 
their data leaving a rms scatter of 9.2 m s _1 (Noyes et al. 1997). More recent data indicate that the orbit 
may have an eccentricity of 0.15 ± 0.03 (Marcy 1998). 

We find the metallicity [Fe/H] = —0.29 ± 0.05 quoted by Gonzalez (1998a) to be only marginally 
consistent with our models. We find consistent models only at the upper end of the quoted error bar. 
Using the observations of Gonzalez (1998a), but with a slightly higher [Fc/H]= — 0.23± 0.06, we calculate a 
mass of 0.89±^M Q , radius of l.3H±° f s R & , and age of 14. lt^ Gyr. We calculate M ce ~ 0.033^^09 and 
t c ~ 21. 51^8- We find general agreement between our computed logg = 4.131q;qq and the observed values 
(See Table 12). 

The activity-age relation predicts an age of 6 Gyr (Noyes et al. 1997), significantly younger than in our 
models. While our old age is worrisome, it is not extremely sensitive to any of the input parameters. For 
example, insisting that the age is <J lOGyr, would require a 2 — a error in both the temperature and the 
metallicity. No rotation has been observed, but based on the Ca II flux, the rotation period is predicted to 
be ~ 20 d (Noyes et al. 1997). Combining this period and our radius, we compute v cq ~ 3.4 km s _1 . This 
is consistent with the observed values of usini ~ 1.5 km s , but an accurate estimate of sini cannot be 
obtained given the large uncertainty in the rotation period. 

4. Implications for Planet Formation and Evolution 

4.1. Tidal Dissipation 

Our models do not change the main conclusions of Rasio et al. (1996) concerning the importance of 
orbital decay driven by tidal dissipation in the 51 Peg system, but they do alter some of the quantitative 
estimates (See Fig. 1). The planets around r Boo and v And have longer timescales for orbital decay, since, 
as our models show, these stars have at most a very shallow convective envelope. The larger semi-major 
axes of the planets around p 1 Cnc and p CrB increase their timescale for orbital decay. Thus, our models 
show that all the known extrasolar planets have stable orbits, in the sense that the orbital decay timescale 
is long compared to the main-sequence lifetime of their star. 

Since orbital circularization is dominated by tidal dissipation in the planet, the only relevant stellar 
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parameter for circularization is the stellar mass, as it affects our determination of the mass of the planet 
and the semi-major axis from observations. Parameters for the planet, such as the radius and the tidal 
dissipation factor, Q, are less certain, but we can estimate them using Jupiter as a guide (Rasio et al. 1996; 
Lubow, Livio, & Tout 1997). Tidal dissipation in the planet could circularize an eccentric orbit for the 51 
Peg, r Boo, and v And systems, and perhaps the p 1 Cnc system as well. Tidal dissipation in the planet 
could not have circularized an eccentric orbit in the case of p CrB. 

Tidal dissipation in the star also tends to spin up the star towards synchronization with the orbital 
motion of the planet, since radial velocity surveys can only determine the vsini, some of the detected 
systems may actually contain low-mass stellar companions in orbits that happen to be viewed nearly 
face-on (sinz <C 1). However, the probability of finding a 1M@ companion with msini as low as 1 Mj or 
10 Mj is only one in ~ 2 • 10 6 or ~ 2 • 10 4 , respectively. Since the SFSU planet search has already found 
6 systems by monitoring a target list of only 120 stars, the spectroscopic binary explanation is extremely 
unlikely, but cannot be strictly ruled out. In principle, astrometry could be used to constrain the maximum 
companion mass. However, this is extremely difficult for short-period systems, since the position wobble is 
so small. For example, assuming a low mass companion to 51 Peg, the amplitude of the induced wobble is 
only 3 • 10 _3 /sini mas. Sophisticated space-based interferometers will be necessary to measure such small 
effects. Absence of X-ray activity can also be used to confirm that a companion is substellar, since close 
spectroscopic binaries are usually observed to have significant X-ray luminosities. However, a young stellar 
age or rapid rotation also correlate with strength of X-ray emission (see e.g., Pravdo et al. 1996). 

Here we focus on the tidal constraints, which can place an upper limit on the mass of companions in 
short-period systems. For the planet to survive, the orbital decay timcscalc must be large compared to the 
age of the system. If the companion has not yet spun up the stellar rotation period to match its orbital 
period, then another constraint can be imposed. The maximum companion mass, as determined from the 
timescales for either orbital decay or spin-up, can be calculated from the stellar parameters. Following 
standard tidal dissipation theory (e.g., Zahn 1977; Zahn & Bouchet 1989; Rasio et al. 1996), we calculate 
the timescale for a planet to spin up the whole starQ as 

* k*-M*R* / n \ 

1 /2 

where = (l*/M*R*\ is the dimensionless gyration radius, is the mass of the star, i?* is the stellar 
radius, m is the mass of the planet, a is the semi-major axis, and r a is the timescale for orbital decay, 
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Here q = m/M* is the mass ratio, / is a numerical factor of order unity, and t c is the eddy turnover 
timescale. Following Rasio et al. 1996, we estimate r c by 

(3) 
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where is the stellar luminosity, M cc is the mass of the convective envelope, and R co is the radius at the 
base of the convective envelope. When the tidal pumping period (half the orbital period) is small compared 



3 Note the obvious typographical error in Eq. 2 of Rasio et al. (1996). Here we use standard tidal dissipation theory based on 
the weak friction approximation with eddy viscosity as the dominant dissipation mechanism in stars with convective envelopes 
(Zahn 1977; Rasio et al. 1996). Although the theory remains controversial (see, e.g., Goodman & Dickson 1998; Tassoul & 
Tassoul 1997), there is empirical support for some of its main predictions (Verbunt & Phinney 1995; Zahn 1992). 



to the eddy turnover time in the convective envelope, we expect that the efficiency of the tidal dissipation 
will be reduced. Thus, for fast tidal pumping, we reduce the factor / according to 
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Since the exact form of this correction is rather uncertain (see Goodman & Oh f997 and references therein), 
we compute tidal decay and spin- up timescales using both / = 1 and / as given by Eq. 4 with /' = 1 . 

An even stronger constraint can be imposed if the spin up of the convective envelope is assumed to 
occur independently from the interior. In the limit of a thin convective shell, the spin-up timescale becomes 



In reality, the true spin- up timescale is likely to be somewhere between and t* u . 
We now consider the importance of tidal effects in each system individually. 

51 Peg. It has been suggested that the 51 Peg system could contain a stellar-mass companion, but 
that the system is being viewed nearly face-on (Kubat, Holmgren, & Rentzsch-Holm 1998; Imbert & Prevot 
1998). Astrometric measurements with Hipparcos impose an upper limit of ~ 500Mj ~ O.5M (Perryman 
ct al. 1997). Diffraction limited (~ 0.05 arcsec) infrared speckle imaging with the Keck telescope searched 
mainly for M dwarf companions and put limits on any such companions. These detection limits give a range 
of maximum orbital separations from 0.75 AU for M3 dwarfs to 23 AU for M7 dwarfs (Marcy et al. 1997). 
The Palomar Testbed Interferometer has also searched the 51 Peg system and determined a maximum 
companion mass of 0.22M Q (Boden et al. 1998). Low levels of X-ray emission also indicate that the 51 Peg 
system is unlikely have a stellar-mass companion (Pravdo et al. 1996). 

We can put another constraint on the maximum companion mass based on the slow rotation rate of 
the star implying t su » tms- We calculate r* u sin 2 z ~ 8.4 • 10 13 yr. Comparing this to the best model 
age of 51 Peg, we get a maximum companion mass of AQMj. If we set / = 1, ignoring the reduction in 
efficiency from the rapid pumping period, then we obtain r* u sin 2 i ~ 1.1 • 10 12 yr, which yields a maximum 
mass of 5.5Mj. We also calculate r s c ° sin 2 i ~ 1.2 • 10 13 yr and t s c „ sin 2 i ~ 1.6 • 10 11 yr, with corresponding 
maximum companion masses of 19Mj and 2.1Mj, respectively (See Fig. 2). 

r Boo. The most interesting result from our models of r Boo is the likely absence of a convective 
envelope. We find convective envelopes only in models with T c ff ^ 6350 K or [Fe/H]^ 0.07. Even in these 
models, the convective envelope is extremely thin. The spectral line bisectors seen by Hatzes et al. (1997) 
and others display a slight curvature, consistent with a very thin convective layer, which our models might 
not capture and which would not be significant for tidal dissipation. Even if we include our models with the 
most significant convective envelopes (T c g — 6200 K, [Fc/H]=0.16), we can place an upper limit on the size 
of the convective zone of M ce <; 0.003M Q (See Table 3). 

While this result does not affect the possibility of tidal dissipation in the planet circularizing the orbit 
as discussed in Rasio et al. (1996) and Rasio & Ford (1996), it does have significant implications for the 
spin history of r Boo. Marcy et al. (1997) and Drake et al. (1998) suggest that the planet is likely to have 
tidally spun up the star. Our models indicate that r Boo does not have a convective envelope and thus 
would not have been spun up by a close companion. In fact, a rotation period of 3.3 d is entirely consistent 
the observed normal rotation rates of young F7 stars with fully radiative envelopes (Gray & Nagor 1985). 
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Thus, the approximate equality of the observed rotation period and the orbital period is most likely a 
coincidence and not the result of tidal spin-up. 

If we insist that tidal dissipation synchronized the stellar rotation period, then it must have occurred 
much earlier in t Boo's life, when the star had a large convective envelope. Our models indicate that r 
Boo's convective envelope starts to shrink rapidly near ~ 2 • 10 6 yr and is reduced to <; 4 • 1CP 3 M after 
~ 1.5 • 10 7 yr. For a short time, ~ 2 — 3 • 10 7 yr, the synchronization timescale is ~ 2 • 10 6 yr, providing a 
window of opportunity for tidal synchronization. By ~ 3 • 10 7 yr, the convective envelope has shrunk to 
less than 1O _4 M . Thus, tidal synchronization would require that the planet be already formed and in a 
short-period orbit before this time. Indeed, it has been suggested that the planet may have migrated to 
its present orbit and stopped at an orbital period synchronized with the stellar rotation, because of the 
tidal force exerted on the orbit (Lin et al. 1996; Trilling et al. 1998). Since the system is quite young, it is 
possible that the star has not had enough time to spin down significantly after the disappearance of the 
convective envelope. 

v And. There is no measured rotation period for v And. Although most of our models of v And do 
have a small convective envelope, some of the models we constructed for T c ff = 6350-R" do not. If we adopt 
our best model, we find that the constraint r* u < tms would impose a limit of 43Mj on the companion 
mass, while ignoring the reduction in efficiency from the rapid pumping period would impose a limit of 
15Mj. The constraint that r s c ° < tms would impose limits of 5.5Mj for /' = 1 and 1.9Mj for / = 1. 

The latest radial velocity data indicate that the planet around v And may have a significantly eccentric 
orbit, e ~ 0.1 (Marcy 1998). If this eccentricity is confirmed, then the non-circular orbit, combined with the 
survival of the planet, could place constraints on the planet's structure. Following Rasio et al. (1996) we 
calculate the timescales for synchronization of the planet's rotation and for orbital circularization. Assuming 
that the planet is identical to Jupiter, we find that the rotation of the planet should be synchronized with 
the orbit and that the circularization timescales is comparable to the age of the system. Since both these 
timescale are very sensitive to the radius of the planet, we can turn the observed non-zero eccentricity into 
an upper limit on the planet's radius. Taking the age to be 2.6 Gyr and still assuming the same dissipation 
factor, Q, as Jupiter, we find R p <J 1.4i?j. Unfortunately, the uncertainty in the age of the system, and 
even more in Q, prevents this from being a stringent constraint at this time. 

p 1 Cnc: We were not able to construct any models of p 1 Cnc with a temperature as low as that 
observed by Gonzalez (1998ab). One possible explanation is that the p 1 Cnc system might be an unresolved 
stellar binary viewed nearly face-on (Gonzalez 1998ab). 

Since the star has an observed rotation period of 42 d, which is longer than the orbital period, we 
can calculate the maximum companion mass which would not have synchronized the stellar rotation 
period. We calculate r* u sin 2 z ~ 1.4 • 10 15 yr. Comparing this to the age of p 1 Cnc, we get a maximum 
companion mass of ~ 0.33M Q . If we set f — 1, ignoring the reduction in efficiency from the rapid pumping 
period, then we obtain r*„ sin 2 i ~ 1.1 ■ 10 14 yr, which yields a maximum mass of 91Mj. We also calculate 
Tg® sin 2 i ~ 3.5 • 10 14 yr and r s c ° sin 2 i ~ 2.5 • 10 13 yr, which give maximum companion masses of ~ 0.1 7M© 
and 46Mj, respectively (See Fig. 4). While this excludes the possibility of most stellar mass companions, we 
can not completely exclude the possibility of a nearly face-on stellar companion of mass O.3M (Gonzalez 
1998ab). Additionally, since the companion is very close to p 1 Cnc, it may be contaminating the spectra 
and distorting the observed parameters sufficiently to affect the maximum mass we calculated here. 

p CrB. While tidal dissipation in the planets could circularize the orbits in 51 Peg, r Boo, v And, 
and possibly also p 1 Cnc, the planet around p CrB has too large a semi-major axis for tidal dissipation to 
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be significant anywhere. Indeed, the orbit is observed to have an eccentricity of 0.15 ± 0.03 (Marcy 1998). 
Depending on the orbital parameters of the secondary star, secular perturbations may be able to produce 
the observed eccentricity in the orbit of the lower mass companion. Alternatively this eccentricity may 
have been produced by a dynamical interaction with another planet. Alternatively, an interaction with a 
circumstellar disk could induce an eccentricity, but only if the planet is significantly more massive than its 
minimum mass (Artymowicz 1992). 

4.2. Metallicity Enhancements 

4-2.1. Observational Evidence 

Observed metallicities are listed in Table 2. Of the five short period planets, all but p CrB have 
[Fe/H]> 0.15. These five stars have a mean metallicity of +0.18. If we consider the four stars with planets 
in the shortest period orbits, then the average metallicity increases to +0.29. Additionally, six of the seven 
stars known to have companions with m sin i <, 5Mj are metal rich relative to the sun and they have a 
mean metallicity of +0.14. These observations have led to speculations that there may be a relationship 
between stars with higher metallicities and stars with planets. 

When analyzing these metallicities, one must be careful to consider the survey population from which 
the stars were selected. As an example, consider the SFSU planet search, which has monitored 120 stars. 
Of those, 67 stars are also in the Cayrel de Strobel et al. (1997) catalogue. Looking at the intersection of 
these two lists, we find that only 8% have an average [Fe/H] > 0.15.0 Since all the stars known to harbor 
planets are F and G stars and it is more difficult to obtain accurate metallicities for K and M stars, it 
may be better to restrict our attention to the 61 F and G stars in the SFSU planet search. Of those, 47 
appear in the Cayrel de Strobel et al. (1997) catalogue, and only about 5% have an average [Fe/H] > 0.15. 
Both populations have a mean [Fe/H]= —0.09 with a standard deviation of 0.22. Additional evidence for 
the a relationship between high metallicities and extrasolar planets comes from a comparison with the 
recent, nearly volume-limited metallicity survey of nearby solar-type stars by Favata et al. (1997). They 
find a mean [Fe/H] of —0.07 ± 0.26 for their full sample and —0.12 ± 0.27 for the portion of the sample 
with T c ff > 5000 K. After they attempt to correct for various biases, they find a peak at [Fe/H] ~ —0.23 
(Favata, Micela, & Sciortino 1997). Extrasolar planets in wide binary systems for which the metallicity 
of each star can be accurately measured could provide useful information about the relationship between 
metallicity enhancements and planet formation. In the cases of r Boo, p 1 Cnc, v And, and p CrB accurate 
determinations of the metallicity of the binary companion are not available. In the 16 Cygni system, 
spectroscopic analyses have revealed very similar metallicities for the two stars, but slight differences in 
the lithium and beryllium abundances (Gonzalez 1998; Lopez & Taoro 1998). Since both these elements 
are involved in nuclear reactions, complications such as rotation can significantly affect the observed 
abundances. Thus, we cannot conclude that there is a significant difference without more careful study. 
As more extrasolar planets are discovered in wide binary systems, these comparisons may become very 
instructive. 

Several factors make some stars better candidates for planet detections than others and introduce 
additional observational biases. For example, high chromospheric activity increases the scatter in the radial 



4 Cayrel de Strobel et al. (1997) catalogues observations of [Fe/H] from the literature. When the catalogue lists multiple 
observations for a single star, we take the unweighted average. 
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velocity measurements, decreasing the sensitivity of the planet search. Some planet searches intentionally 
exclude such stars, while others observe them in spite of the reduced detection efficiency. In either case, 
this introduces complicated selection effects. While the selection effects arc not fully understood, it seems 
that high stellar mctallicities may be correlated with the existence of planets, or at least with short period 
planets. 

4-2.2. Mechanisms to Accrete High-Z Material 

Several mechanisms that have been discussed in the context of extrasolar planet formation may lead to 
the accretion of high-Z material onto the parent star. While such theories might also explain the apparent 
correlation of short period planets with higher metallicities, the plausibility of such an explanation depends 
on many factors. 

One factor is the type of material accreted onto the star. Accreting material from a gaseous disk will 
have little effect, since this material has essentially the same metallicity as material which formed the 
star. However, accreting a large number of terrestrial planets, asteroids, or planetesimals could produce a 
noticeable effect. Alternatively, accreting gaseous giant planets like Jupiter could also have a significant 
effect since these planets are thought to contain up to ~ 10M e of rocky material in their core. 

The size of the stellar convective envelope at the time of accretion is also of critical importance, 
since the accreted material is rapidly diluted over the entire stellar convective region. In particular, if 
the accretion of high-Z material occurred while the star was still on the PMS and had a large convective 
envelope, then there would be little effect on the stars' observed metallicities. To be efficient in increasing 
the surface metallicity, the accretion must take place sufficiently late in the stellar evolution when the outer 
convective envelope is shallow (See Fig. 5). 

Finally, the particular type of accretion process which adds high-Z material to the star could also be 
important. If small amounts of high-Z material are gradually deposited from nearly circular orbits, then it 
might remain in the outer convective envelope. However, if a massive object entered the star from a highly 
eccentric orbit, it might penetrate more deeply into the star, thereby diluting the high-Z material over a 
larger mass fraction. Furthermore, if a rocky body is added to the star from a nearly circular orbit, it still 
might sink to a significant depth before being disrupted. Even if the high-Z material is not immediately 
mixed deep into the stellar interior, a layer of high-Z material would have a greater mean molecular weight, 
which could drive a thermal instability allowing it to penetrate deeper into the star. It is beyond the scope 
of this paper to take all of these process into account. Instead, we simply obtain an upper bound on the 
effects of adding high-Z material to the star by noting that any material added must have been diluted at 
least across the present convective zone. 

Our models show that, at present, 51 Peg, p 1 Cnc, and p CrB all have significant convective envelopes, 
while v And and r Boo have either very thin convective envelopes or none at all (See Table 3). However, 
our best models of v And and r Boo indicate that they did have convective envelopes at the ZAMS (See 
Table 3). Since the size of the current convective envelope imposes an upper limit on the effects of chemical 
pollution from accretion independent of when the accretion occurred, we will only consider 51 Peg, p 1 Cnc, 
and p CrB. We discuss the effects of adding various high-Z materials to these stars, but with an original 
metallicity of [Fe/H]~ 0. This metallicity is already higher than most nearby solar-type stars, as discussed 
above. 
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Gas Giants. Lin ct al. (1996) proposed a model in which several multiple Jupiter-mass planets 
migrate inwards as a result of dissipation in a massive disk. Today, we are only able to observe the last of 
the planets which were fortunate enough to barely avoid spiraling into their parent star. We consider the 
effect of adding Jupiter-like giant planets with 10M ffi cores of roughly chondritic composition (Anders & 
Grevesse 1989). Starting with a solar metallicity star and the presently observed convective envelope, it 
would take ~ 5 and ~ 15 such Jupiter-like giant planets to raise the stellar metallicities to those observed 
in 51 Peg and p 1 Cnc, respectively. However, in this model, the accretion takes place while the star still has 
a gaseous protostellar disk (t <, 10 7 yr). At this stage the star would still be on the PMS where it would 
have a much larger convective envelope. Therefore, much more mass would be necessary to raise their 
metallicities to those observed today (See Fig. 5). 

In a similar model, a single massive giant planet (<~ 4Mj) migrates towards the parent star and then 
undergoes Roche lobe overflow, stabilizing its orbit against spiraling further inward (Trilling ct al. 1998). 
Again, gas giant material would be accreted onto the star, although this time the rocky core would not be 
accreted. Since Jupiter's atmosphere is only slightly enriched in heavy elements, this is not likely to alter a 
star's observed metallicity significantly. Additionally, if the migration is due to a viscous disk, the star will 
still be highly convective, rendering the accreted material inefficient in raising the surface metallicity. 

Gaseous Disks: In both of the above models, a portion of the gas disk can also be accreted onto the 
star. The disk may have a metallicity slightly larger than the star, but this is not expected to significantly 
alter the observed stellar metallicity Since both the disk mass and the disk metallicity will vary with time 
and neither are well known, a quantitative estimate for this mechanism would be difficult. 

In contrast, adding rocky material (terrestrial planets, asteroids, or planetesimals) to the star is much 
more efficient at increasing the observed values of [Fe/H], as they have a much higher metallicity than 
gaseous giant planets. Given the presently observed convective envelopes and an initially solar composition, 
it would require only ~ 40M & and ~ IOOM0 of chondritic material (Anders & Grevesse 1989) to reach the 
observed metallicities of 51 Peg and p 1 Cnc, respectively. However, standard planet formation scenarios do 
not predict much more than ~ 20M ffi in terrestrial planets. 

Planetesimals: Accretion of a large amount of solid material in the form of asteroids or planetesimals 
arises naturally in the model proposed by Murray et al. (1998). In this model giant planets can migrate 
inwards while scattering planetesimals into the parent star. 

Since the early migration is rapid, the planet will not clear out the planetesimal disk until its migration 
slows. Thus we neglect any planetesimals that are scattered into the star until the planet reaches its final 
orbit. The mass accreted on the star is given by 



is the fraction of planetesimals scattered onto the star by the planet at distance a, m is the mass of the 
planet, and a is a parameter ~ 0.5 — 1.0, which must be weighted across multiple resonances (Murray et 
al. 1998; Hansen 1998). Note that /(a) only accounts for planetesimals which become planet crossing before 
colliding with the star, i.e., planetesimals which collide with the star without ever becoming planet-crossing 
are neglected. 



M, 



acc 




(6) 



where 




(7) 



Using this prescription, we can estimate the mass of planetesimals which would be scattered into the 
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star for 51 Peg, p 1 Cnc and p CrB. We find 0.4Mj (130M e ) for 51 Peg and 0.6Mj (180M e ) for both p 1 
Cnc, and p CrB. Since all this mass is in asteroids, this could lead to a significant and observable increase 
in the metallicity. Starting from a star with solar metallicity and the present convective envelope, adding 
130M ffi of asteroids to 51 Peg or 180M e of asteroids to p 1 Cnc increases the observed [Fe/H] to 0.48 and 
0.39 dex, respectively, assuming (our best model for the two stars). However, the requirement that the 
planetesimal disk still be present, implies that the stars are still young with convective envelopes larger than 
at present. The time it takes for a planet to migrate to a 0.05 AU orbit in the Murray et al. (1998) scenario 
is <, 3 • 10 7 yr, but this depends on the model of the planetesimal disk. Once the planet stops migrating, 
the asteroids inside its orbit are quickly cleared out. As can be seen from Fig. 5, there is a window of 
opportunity at 2 — 3 • 10 yr which is the right time to produce the observed metallicities. However, since 
the convective envelope disappears rapidly, this requires a very specific timing. Different models of the disk 
or the inclusion of additional resonances could destroy any careful tuning of parameters in a particular 
model adjusted to produce significant metallicities. 

Terrestrial Planets: Rasio & Ford (1996) and Weidenschilling & Mazari (1996) proposed that 51 
Peg- like planets could be produced by gravitational scattering off a second massive planet. If a significant 
number of planets are to wind up very close to their parent stars, a significant number should also have 
collided with their parent stars. However, it is unlikely that two planets in the same system would be 
scattered so close to their parent star. Thus, one would not expect to see increased metallicities in stars 
with 51 Peg-like planets due to the accretion of giant planets. Similarly, in the model of Kiseleva & Eggleton 
(1997), secular perturbations of a binary companion induce large eccentricities which are later circularized. 
In either this model or the previous model, a giant planet acquires a very high eccentricity which would 
disrupt inner terrestrial planets, possibly making them collide with the parent star. These mechanisms 
have the advantage that they can deposit material onto the star after it has reached the main sequence 
and its convective envelope has diminished to near its present size. It would take a large amount of mass 
in terrestrial planets and asteroids for these mechanisms to increase the metallicity from solar to what is 
observed in 51 Peg or p 1 Cnc (~ 40M e and ~ 100M ffi , respectively). One would still expect the planet to 
scatter asteroids into the star, as in the Murray et al. (1998) scenario. However, here a particular disk mass 
or timing is not required. 

4.2.3. Polluted Stellar Models 

The possibility of a star having a surface metallicity significantly higher than its interior metallicity 
arises naturally in the context of planet formation. Thus, when constructing models based on observational 
data, we should realize that the chemical composition of the interior, which contains most of the stellar 
mass, could deviate significantly from the composition observed on the stellar surface. For this reason, it is 
important to consider stellar models with interior metallicities somewhat lower than those observed. 

We have constructed a few models to explore the effects of adding high-Z material to the surface of a 
star. We pick an initial metallicity for the entire star and evolve a stellar model slightly past the ZAMS. 
Then, we increase the metallicity in the outer convective envelope and continue to evolve the star on the 
main sequence. 

First, we consider the effect of adding a Jupiter with 10M ffi of high-Z material (in its rocky core, using 
the procedure described above) to the ZAMS Sun. We find that the effective temperature, luminosity, and 
radius increase by 0.55%, 2.6%, and 0.19%, respectively, while the mass of the convective envelope decreases 
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by approximately 5.5%. Similar models have been investigated in the context of the Solar neutrino problem 
and have been found to reduce the agreement between theoretical and observational determinations of solar 
p-mode frequencies (Jeffery, Mailey, & Chambers 1997; Christensen-Dalsgaard & Gough 1998; Bahcall 
1989). 

Next, we consider the consequences of a reduced interior metallicity for our models. Since we use 
observations of the luminosity, effective temperature, and surface metallicity, the radius will not be affected, 
but the mass, age, and size of the convective envelope will be. Generally, decreasing the interior metallicity 
from the observed surface metallicity causes the mass to decrease, but the age and the size of the convective 
envelope both increase. For example, if we take the observed parameters for 51 Peg, but set the interior 
metallicity to solar, then the we obtain in our new best model a mass of 0.94M©, an age of 10.8 Gyr, and a 
convective envelope of mass 0.026M Q (see Table 3 for comparison). 

Finally, we have explored the possibility that a reduced interior metallicity could explain the difficulties 
encountered in constructing models that match the observed parameters for p 1 Cnc and r Boo. In general, 
lowering the interior metallicity leads to a higher surface temperature in our models. Therefore, in the 
case of p 1 Cnc, this would make the discrepancy worse. However, this could bring our models of r Boo 
closer to agreement with the observed effective temperature. If we demand that our models of r Boo have 
T c ff = 6550 K, as indicated by the recent observations of Gonzalez (1998a), then the difference between the 
interior and surface metallicities must be ~ 0.35 dex. Although possibly a coincidence, this is comparable 
to t Boo's observed [Fe/H]=0.34 ± 0.09. However, since we were able to construct models with a single 
metallicity for T e g = 6500 K, the uncertainty in the measurement of the effective temperature (~ 100 K) 
prevents us from placing a firm constraint on the interior metallicity. If we use the same parameters for r 
Boo as discussed in §3.2, but with a solar metallicity interior, we find a mass of 1.29M©, an age of 0.4 Gyr, 
and no significant outer convective envelope. 

We find that a reduced interior metallicity could resolve the difficulty in constructing models that 
match the observed luminosity, temperature, and metallicity for r Boo. In our calculations we have 
instantaneously distributed the accreted material across the convective envelope of the star and then 
evolved the star in this dual composition state. If the material is deposited on a dynamical timescale, then 
the star will temporarily be out of thermal equilibrium. A future calculation could evolve the star on a 
thermal timescale to follow the readjustment to a new thermal equilibrium. In addition, if the convective 
envelope has a higher metallicity than the interior, then the possibility of triggering a thermal instability 
that could drive high metallicity material into the stellar interior should be examined. 

5. Summary 

We have computed stellar models for nearby stars with close planetary-mass companions. Tables 5, 
7, 9, 11, and 13 show the properties of models constructed for different values of the observed parameters 
within measurement uncertainties. These models were used to calculate the stellar parameters and error 
bars summarized in Table 3. We will make a larger grid of stellar models available electronically to help 
interpret future stellar observations. 

Using conventional tidal dissipation theory, we have studied tidal dissipation effects in these systems, 
based on the results of our models. We find that the orbital decay timescale is longer than the age of the 
star in all systems (See Figs. 1 and 3). The timescale for orbital circularization is shorter than the stellar 
ages in 51 Peg, r Boo, and v And, comparable to the stellar age in p 1 Cnc, but longer than the stellar 
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age in p CrB. We have also calculated maximum companion masses based on the lack of synchronization 
between the stellar rotation and the companion orbital period for 51 Peg, v And and p 1 Cnc (See Figs. 2 
and 4). For the most optimistic assumptions for the tidal coupling (in particular assuming that the tidal 
torque is exerted only on the outer convective envelope), we find m <, 2.1Mj, 1.9 Mj and 46 Mj for 51 
Peg, v And, and p 1 Cnc, respectively. For r Boo the similarity between the stellar rotation period and the 
orbital period is likely a coincidence and not the result of tidal spin-up. If the significant eccentricity of the 
v And system is confirmed, then we can place an upper limit on the radius of the planet, R p <J 1.4i?,j up . 

We have examined the observational evidence for a correlation between the presence of a planet and a 
high stellar metallicity. We have evaluated the effects of accreting high-metallicity material onto a star in 
the context of various proposed mechanisms for producing short period systems. The accretion of gas giants 
or of material from a gaseous disk causes only minor metallicity enhancements if the stellar convective 
envelope is still large at the time of the accretion. Accretion of rocky material, such as planetesimals, 
asteroids, and terrestrial planets, could cause significant increases in observed metallicities if enough mass 
is available. An increased surface metallicity can also have a significant effect on stellar models and on the 
interpretation of observational data. 
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Fig. 1. — This plot of planet mass versus orbital period indicates the regions allowed after imposing tidal 
stability constraints on 51 Pegasi. The lines labeled r a indicate where the timescale for orbital decay of the 
planet would be less than the age of the star, according to standard tidal dissipation theory. The lines labeled 
T' a differ in that no reduction of the orbital decay timescale was applied to account for the long convective 
turnover timescale (Eq. 4). For both groups of lines labeled T a and T' a , the thicker line corresponds to our 
best model and the lines on either side reflect the uncertainty in the stellar age from our models (see Tabic 
15). The lines labeled K/ smi indicate the typical detection limits of current radial velocity surveys. The 
line labeled Rl/Rj indicates where a planet with mass 0.4Mj, and radius ~ lRj ~ O.li?©, would overflow 
its Roche lobe. The large dot indicates the minimum mass of the planet. This plot was constructed using 
our best model parameters: M* = 1.05M Q , i?» = 1.16.R©, t c = 18.6 d, M cc = 0.023M Q , and t ms = 7.6±*-° 
Gyr. 

Fig. 2. — This plot shows several time scales as a function of the mass of the companion orbiting 51 Peg. 
Here r a is the orbital decay timescale, r* u is the spin-up timescale for the whole star, and is the spin-up 
timescale for the convective envelope alone. The dashed diagonal lines correspond to setting / = 1, while the 
solid diagonal lines correspond to calculating / from Eq. 4 with /' = 1. The solid horizontal line indicates 
the age of the system, tms, as determined from our best model, while the dashed horizontal lines indicate 
the uncertainty in the age. The intersections of the age with the spin-up timescales gives the maximum 
companion masses quoted in the text. 

Fig. 3. — This plot of planet mass versus orbital period indicates the regions allowed by tidal stability 
constraints on r Boo. Here we have taken our model with the largest convective envelope. The model has 
M* = 1.3OM , R* = 1.51.R©, t c = 7.93 d, M cc = O.OO35M , and an age of 2.9 Gyr. This model is for 
T off = 6200 K, [Fc/H]= 0.16, and L* = 3.006X©. Other conventions are as in Fig. 1. 

Fig. 4. — This plot shows several time scales as a function of the mass of the companion orbiting p 1 Cnc. 
The labels are the same as in Fig. 2. The intersection of the age with the timescales for tidal spin- up lead 
to the maximum companion masses listed in the text. 

Fig. 5. — The solid line (left axis) shows the evolution of the Sun's convective envelope. The dashed lines 
(right axis) indicate the surface metallicity which would result from the instantaneous accretion of rocky 
material onto the star at each time in the Sun's past, assuming that the accreted material is mixed across 
the convective envelope. Starting from the bottom, the dashed lines correspond to accreting 5, 10, 25, 50, 
100, 500, and 1000 M e . We see that producing the high surface mctallicities seen in 51 Peg (0.21 dex) and 
p 1 Cnc (0.29 dex) is possible with the accretion of a large amount of rocky material after ~ 10 7 yr. 
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Table 1. Dynamical Properties of Extrasolar Planets a 



Star 


HD 


M. 


Porb 


e 


a 


m sin i 


Ref 






(M ) 


(days) 




(AU) 


(Mj) 




r Boo 


120136 


1.37 ± 0.09 


3.3128 ± 0.0002 


0.018 ± 0.016 


0.0413 ± 0.0010 


4.74 ± 0.10 


[1] 


51 Peg 


217014 


1.05 ± 0.09 


4.2311 ± 0.0005 


0.012 ± 0.010 


0.0520 ± 0.0015 


0.464 ± 0.036 


[2] 


v And 


9826 


1.34 ± 0.12 


4.611 ± 0.005 


0.109 ± 0.040 


0.0598 ± 0.018 


0.739 ± 0.075 


[1] 


p 1 Cnc 


75732 


0.95 ± 0.10 


14.648 ± 0.0009 


0.051 ± 0.013 


0.1152 ± 0.0038 


0.874 ± 0.079 


[1] 


p CrB 


143761 


0.89 ± 0.05 


39.645 ± 0.088 


0.15 ± 0.03 


0.2192 ± 0.0042 


1.064 ± 0.059 


[3,7] 


HD 114762 


114762 


0.75 ± 0.15 


83.90 ± 0.08 


0.35 ± 0.05 


0.34 ± 0.04 


10.5 ± 1.9 


[4] 


70 Vir 


117176 


1 01+ 005 

-0.02 


116.67 ± 0.01 


0.4 ± 0.01 


0.47 ± 0.01 


7.3 ± 0.3 


[5] 


16 Cyg B 


186427 


0.96 ± 0.05 


800.8 ± 11.7 


0.634 ± 0.082 


1.66 ± 0.05 


2.0 ± 0.3 


[6] 


47 UMa 


95128 


1.01 ± 0.05 


1090 ± 15 


0.03 ± 0.06 


2.08 ± 0.06 


2.3 ± 0.2 


[5] 



a The orbital parameters (orbital period P rb, eccentricity e, semimajor axis a, and minimum 
mass msin?) are from: [1] Butler & Marcy 1996; [2] Marcy et al. 1997; [3] Noyes et al. 1997; 
[4] Mazeh et al. 1996; [5] Marcy & Butler 1996; [6] Cochran et al. 1997a; [7] Marcy 1998. In 
order to determine a and msini, we have assumed the stellar mass given in column 3. The first 
five stellar masses are based on the results of this paper and are discussed later. The remaining 
stellar masses are based on preliminary calculations which will be presented in a subsequent 
paper. 
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Table 2. Observed Parameters of Stars with Extrasolar Planets b 



Star 


MK 


M v 


Teff (K) 


[Fc/H] 


log g (cgs) 


BC (V) 


L,/L @ 


r Boo 


F7V 


3.535 ± 0.024 


6550 ± 100 


+0.34 ± 0.09 


4.50 ± 0.15 


+0.02 


3.006 ± 0.213 


51 Peg 


G5V 


4.518 ± 0.025 


5750 ± 75 


+0.21 ± 0.06 


4.40 ± 0.10 


-0.07 


1.321 ± 0.094 


v And 


F8V 


3.453 ± 0.021 


6250 ± 100 


+0.17 ± 0.08 


4.30 ± 0.10 


-0.01 


3.333 ± 0.225 


p 1 Cnc 


G8V 


5.470 ± 0.023 


5250 ± 70 


+0.45 ± 0.03 


4.40 ± 0.15 


-0.18 


0.608 ± 0.042 


p CrB 


G2V 


4.184 ± 0.027 


5750 ± 75 


-0.29 ± 0.06 


4.10 ± 0.05 


-0.07 


1.796 ± 0.132 


HD 114762 


F9V 


4.259 ± 0.123 


5950 ± 75 


-0.60 ± 0.06 


4.45 ± 0.05 


-0.07 


1.676 ± 0.290 


70 Vir 


G5V 


3.680 ± 0.029 


5500 ± 75 


-0.03 ± 0.06 


3.90 ± 0.10 


-0.07 


2.856 ± 0.076 


16 Cyg B 


G5V 


4.597 ± 0.024 


5700 ± 75 


+0.06 ± 0.06 


4.35 ± 0.05 


-0.07 


1.228 ± 0.087 


47 UMa 


GOV 


4.288 ± 0.020 


5800 ± 75 


+0.01 ± 0.06 


4.25 ± 0.05 


-0.05 


1.603 ± 0.106 



b Spectral type from Hipparcos Catalogue; My derived from Hipparcos Catalogue (ESA, 
1997) photometric and parallax data; T e g, [Fe/H], and log g from Gonzalez 1997ab, 1998ab; 
BC(V) derived from previous parameters using color calibration as in the Yale isochrones 
(Green et al. 1987); calculated from the Hipparcos absolute visual magnitude and 

BC(V). 
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Table 3. Primary Results of Stellar Model Calculations 



Star 


M, 
(Mq) 


R, 

(Re) 


Ago 
(GYr) 




(M ) 


Rcc 

(Re) 


(d) 






Mce.ZAMS 

(M @ ) 


r Boo 


1.37 ± 0.08 


1 4i+oio 

-0.09 


1 2+ 1 ' 2 




^ 0.002 


1.19 ± 0.09 


^ 5.8 


0.10 


£ 0.002 


0.008 


51 Peg 


1 nc+O.OO 


1.16 ± 0.07 


7 6+ 40 

' - u -5.1 





n9 ,+0.007 
u ^°-0.006 


0.82 ± 0.04 


18.6 ± 2.5 


0.12 


0.020 


0.037 


v And 


I Q 4 + 0.07 


1 56+ 011 




0. 


nn9 +0.003 
uuz -0.002 


, 9q + 0.07 


R 8 + 23 

6.8_ 6 8 


0.10 


0.002 


0.006 


p 1 Cnc 


Q5+ 011 


93+ 02 




0. 


04 g+0.004 
"^"-0.006 


0.63 ± 0.03 




0.15 


0.048 


0.069 


p CrB 


89 +0 0s 


, oc+0.09 
1,M -0 .08 


14.1+. 2 ;° 


0. 


„„, + 0.011 
UJO -0.009 


o.88;° ° 5 3 


21.5+- 


0.10 


0.021 


0.029 
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Table 4. Summary of Observational Data for 51 Peg 



T cff (K) 


[Fc/H] (solar) 


log 5 (cgs) 


v sin i (km s 1 ) 


Reference 


5750 ± 75 


0.22 ±0.06 


4.40 ±0.10 


1.4 ±0.3 


Gonzalez 1997a, 1998 


5793 ± 70 


0.20 ±0.07 


4.33 ±0.10 


2.0 ± 1.0 


Fuhrmann et al. 1997 








2.35 ±0.1 


Hatzes et al. 1997 


5775 


0.20 ± 0.07 


4.35 




Tomkin et al. 1997 








2.4 ±0.3 


Francois et al. 1996 


5669 




4.06 




Gratton et al. 1996 




0.172 ± 0.049 






Taylor 1996 c 




0.19 


4.32 


2.8 ±0.5 


Mayor & Queloz 1995 


5730 








Gray 1995 


5773 








Grennon (cited by Mayor & Queloz 1995) 


5724 




4.30 




Valenti (cited by Mayor & Queloz 1995) 


5755 


0.06 


4.18 




Edvardsson 1993b 


5740 


0.05 


3.76 




McWilliam 1990 






4.13 




Xu 1991 


5750 


0.12 




1.7 ±0.8 


Soderblom 1982 


5727 


0.12 


4.27 




Hearnshaw 1972 



c This value is based on a statistical analysis of previous measurements and not any new 
observations. 
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Table 5. Stellar Models for 51 Peg 



M* 


m 


Toff 


L, 


Z 


X 


Y 


Ago 


M ce 


PSH 




L, 


R, 


T {t 


log9 


SY 

Tz 


Solar 


Solar 


K 


Solar 








Gyr 


Solar 


Solar 


Solar 


Solar 


Solar 


K 


cgs 


Solar 



1.092 


0.21 


5750 


1.321 


0.030 


0.687 


0.283 


6.4 


0.02226 


0.0942 


0.835 


1.322 


1. 


161 


5751 


4.35 


0.000 


1.067 


0.21 


5750 


1.321 


0.030 


0.669 


0.301 


7.2 


0.02313 


0.0954 


0.827 


1.321 


1. 


161 


5750 


4.34 


1.500 


1.050 


0.21 


5750 


1.321 


0.030 


0.657 


0.312 


7.6 


0.02371 


0.0962 


0.822 


1.321 


1. 


161 


5750 


4.33 


2.500 


1.033 


0.21 


5750 


1.321 


0.030 


0.646 


0.324 


8.1 


0.02429 


0.0970 


0.816 


1.320 


1. 


161 


5749 


4.32 


3.500 


1.059 


0.21 


5675 


1.321 


0.030 


0.687 


0.283 


8.8 


0.02772 


0.1013 


0.833 


1.321 


1 


192 


5675 


4.31 


0.000 


1.033 


0.21 


5675 


1.321 


0.030 


0.669 


0.301 


9.5 


0.02856 


0.1024 


0.824 


1.321 


1 


192 


5675 


4.30 


1.500 


1.016 


0.21 


5675 


1.321 


0.030 


0.657 


0.312 


10.1 


0.02912 


0.1031 


0.819 


1.321 


1. 


192 


5675 


4.29 


2.500 


0.998 


0.21 


5675 


1.321 


0.030 


0.646 


0.324 


10.6 


0.02968 


0.1039 


0.813 


1.321 


1. 


192 


5675 


4.28 


3.500 


1.124 


0.21 


5825 


1.321 


0.030 


0.687 


0.283 


4.2 


0.01769 


0.0875 


0.836 


1.324 


1. 


132 


5828 


4.38 


0.000 


1.100 


0.21 


5825 


1.321 


0.030 


0.669 


0.301 


4.9 


0.01843 


0.0886 


0.829 


1.324 


1. 


132 


5827 


4.37 


1.500 


1.084 


0.21 


5825 


1.321 


0.030 


0.657 


0.312 


5.3 


0.01893 


0.0893 


0.824 


1.323 


1. 


132 


5827 


4.37 


2.500 


1.068 


0.21 


5825 


1.321 


0.030 


0.646 


0.324 


5.8 


0.01942 


0.0901 


0.819 


1.323 


1.132 


5826 


4.36 


3.500 


1.050 


0.15 


5750 


1.321 


0.027 


0.690 


0.283 


7.8 


0.02328 


0.0960 


0.825 


1.325 


1. 


161 


5755 


4.33 


0.000 


1.032 


0.15 


5750 


1.321 


0.027 


0.679 


0.295 


8.4 


0.02390 


0.0968 


0.819 


1.324 


1. 


161 


5753 


4.32 


1.500 


1.020 


0.15 


5750 


1.321 


0.027 


0.671 


0.302 


8.8 


0.02431 


0.0974 


0.815 


1.323 


1. 


161 


5752 


4.32 


2.500 


1.008 


0.15 


5750 


1.321 


0.027 


0.663 


0.310 


9.2 


0.02472 


0.0980 


0.811 


1.322 


1. 


161 


5751 


4.31 


3.500 


1.099 


0.27 


5750 


1.321 


0.035 


0.658 


0.307 


6.0 


0.02224 


0.0938 


0.835 


1.321 


1. 


161 


5750 


4.35 


1.500 


1.073 


0.27 


5750 


1.321 


0.035 


0.641 


0.324 


6.7 


0.02294 


0.0950 


0.828 


1.321 


1. 


161 


5750 


4.34 


2.500 


1.048 


0.27 


5750 


1.321 


0.035 


0.625 


0.340 


7.5 


0.02364 


0.0962 


0.820 


1.321 


1. 


161 


5749 


4.33 


3.500 


1.092 


0.21 


5750 


1.233 


0.030 


0.687 


0.283 


5.6 


0.02203 


0.0902 


0.812 


1.234 


1. 


122 


5751 


4.38 


0.000 


1.068 


0.21 


5750 


1.233 


0.030 


0.669 


0.301 


6.3 


0.02287 


0.0913 


0.804 


1.233 


1. 


122 


5750 


4.37 


1.500 


1.052 


0.21 


5750 


1.233 


0.030 


0.657 


0.312 


6.8 


0.02343 


0.0920 


0.799 


1.233 


1. 


122 


5750 


4.36 


2.500 


1.036 


0.21 


5750 


1.233 


0.030 


0.646 


0.324 


7.3 


0.02400 


0.0928 


0.794 


1.233 


1. 


122 


5750 


4.35 


3.500 


1.094 


0.21 


5750 


1.415 


0.030 


0.687 


0.283 


7.1 


0.02251 


0.0982 


0.860 


1.415 


1.202 


5750 


4.32 


0.000 


1.068 


0.21 


5750 


1.415 


0.030 


0.669 


0.301 


7.8 


0.02328 


0.0994 


0.852 


1.415 


1.202 


5750 


4.31 


1.500 


1.050 


0.21 


5750 


1.415 


0.030 


0.657 


0.312 


8.2 


0.02379 


0.1002 


0.846 


1.415 


1.202 


5750 


4.30 


2.500 


1.033 


0.21 


5750 


1.415 


0.030 


0.646 


0.324 


8.7 


0.02431 


0.1010 


0.841 


1.415 


1.202 


5750 


4.29 


3.500 
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Table 6. Summary of Observational Data for t Boo 



T e ff (K) [Fe/H] (solar) logg (cgs) vs'mi (km s x ) 



6360 ±80 0.27 ±0.08 4.17 ±0.10 15.6 ±0.7 
6550 ± 100 

6600 ± 100 0.34 ± 0.09 4.50 ±0.15 14.5 ± 0.5 

6405 0.18 

15 ±1 

6358 0.09 4.22 

6390 0.30 3.8 

6460 0.00 4.3 

6300 0.00 4.60 

14.8 ±0.3 

6490 0.21 17 ± 1 

6380 0.14 4.3 

6450 0.28 4.3 

6462 0.28 4.3 



Reference 



Fuhrmann et al. 1998 

Gonzalez 1997b (spectroscopic upon additional analysis) 

Gonzalez 1997a, 1998 (spectroscopic) 

Gonzalez 1997a (photometric) 

Baliunas et al. 1997 

Marsakov et al. 1995 

Boesgaard & Lavery 1986 

Thevenin et al. 1986 

Thevenin & Foy 1983 

Gray 1982 

Soderblom 1982 

Kuroczkin & Wiszniewski 1977 

Perrin 1977 

Spite 1968 
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Table 7. Stellar Models for t Boo 



M* 






L, 


Z 


X 


Y 


Ago 


M ce 


PSH 




L, 


R, 


T a ff 


logs 


SY 

sz 


Solar 


Solar 


K 


Solar 








Gyr 


Solar 


Solar 


Solar 


Solar 


Solar 


K 


cgs 


Solar 



1.403 


0.25 


6400 


2.808 


0.033 


0.684 


0.283 


0.5 


0.00000 


0.0000 


1.154 


2.816 


1.367 


6404 


4.31 


0.000 


1.375 


0.25 


6400 


2.808 


0.033 


0.662 


0.305 


0.7 


0.00000 


0.0000 


1.157 


2.812 


1.367 


6402 


4.30 


1.500 


1.357 


0.25 


6400 


2.808 


0.033 


0.647 


0.320 


0.9 


0.00000 


0.0000 


1.159 


2.809 


1.367 


6401 


4.30 


2.500 


1.338 


0.25 


6400 


2.808 


0.033 


0.632 


0.334 


1.1 


0.00000 


0.0000 


1.161 


2.807 


1.366 


6399 


4.29 


3.500 


1.350 


0.16 


6400 


3.006 


0.027 


0.690 


0.283 


1.5 


0.00000 


0.0000 


1.191 


3.008 


1.414 


6401 


4.27 


0.000 


1.337 


0.16 


6400 


3.006 


0.027 


0.677 


0.296 


1.6 


0.00000 


0.0000 


1.193 


3.007 


1.414 


6400 


4.26 


1.500 


1.329 


0.16 


6400 


3.006 


0.027 


0.669 


0.304 


1.7 


0.00000 


0.0000 


1.194 


3.007 


1.414 


6400 


4.26 


2.500 


1.320 


0.16 


6400 


3.006 


0.027 


0.661 


0.312 


1.7 


0.00000 


0.0000 


1.196 


3.007 


1.414 


6400 


4.26 


3.500 


1.450 


0.34 


6350 


3.006 


0.041 


0.676 


0.283 


0.8 


0.00000 


0.0000 


1.230 


3.015 


1.437 


6353 


4.28 


0.000 


1.416 


0.34 


6350 


3.006 


0.041 


0.642 


0.316 


1.0 


0.00000 


0.0000 


1.224 


3.013 


1.437 


6353 


4.27 


1.500 


1.392 


0.34 


6350 


3.006 


0.041 


0.620 


0.339 


1.1 


0.00024 


0.0112 


1.220 


3.012 


1.437 


6353 


4.27 


2.500 


1.369 


0.34 


6350 


3.006 


0.041 


0.598 


0.361 


1.3 


0.00096 


0.0452 


1.215 


3.011 


1.436 


6352 


4.26 


3.500 


1.383 


0.25 


6300 


3.006 


0.033 


0.684 


0.283 


1.7 


0.00140 


0.0729 


1.241 


3.013 


1.460 


6302 


4.25 


0.000 


1.362 


0.25 


6300 


3.006 


0.033 


0.662 


0.305 


1.8 


0.00163 


0.0750 


1.232 


3.013 


1.460 


6302 


4.24 


1.500 


1.349 


0.25 


6300 


3.006 


0.033 


0.647 


0.320 


1.9 


0.00178 


0.0763 


1.227 


3.013 


1.460 


6303 


4.24 


2.500 


1.335 


0.25 


6300 


3.006 


0.033 


0.632 


0.334 


2.0 


0.00193 


0.0777 


1.221 


3.013 


1.459 


6303 


4.23 


3.500 


1.423 


0.25 


6450 


3.006 


0.033 


0.684 


0.283 


0.5 


0.00000 


0.0000 


1.159 


3.008 


1.392 


6451 


4.30 


0.000 


1.397 


0.25 


6450 


3.006 


0.033 


0.662 


0.305 


0.7 


0.00000 


0.0000 


1.162 


3.008 


1.392 


6451 


4.30 


1.500 


1.380 


0.25 


6450 


3.006 


0.033 


0.647 


0.320 


0.8 


0.00000 


0.0000 


1.164 


3.008 


1.392 


6451 


4.29 


2.500 


1.363 


0.25 


6450 


3.006 


0.033 


0.632 


0.334 


1.0 


0.00000 


0.0000 


1.166 


3.008 


1.392 


6451 


4.28 


3.500 


1.408 


0.25 


6400 


3.006 


0.033 


0.684 


0.283 


0.9 


0.00000 


0.0000 


1.189 


3.005 


1.414 


6400 


4.29 


0.000 


1.384 


0.25 


6400 


3.006 


0.033 


0.662 


0.305 


1.1 


0.00000 


0.0000 


1.192 


3.006 


1.414 


6400 


4.28 


1.500 


1.368 


0.25 


6400 


3.006 


0.033 


0.647 


0.320 


1.2 


0.00000 


0.0000 


1.194 


3.007 


1.414 


6400 


4.27 


2.500 


1.351 


0.25 


6400 


3.006 


0.033 


0.632 


0.334 


1.3 


0.00000 


0.0000 


1.196 


3.008 


1.414 


6401 


4.27 


3.500 


1.419 


0.25 


6400 


3.218 


0.033 


0.684 


0.283 


1.2 


0.00000 


0.0000 


1.224 


3.220 


1.463 


6401 


4.26 


0.000 


1.396 


0.25 


6400 


3.218 


0.033 


0.662 


0.305 


1.3 


0.00000 


0.0000 


1.228 


3.218 


1.463 


6400 


4.25 


1.500 


1.381 


0.25 


6400 


3.218 


0.033 


0.647 


0.320 


1.4 


0.00000 


0.0000 


1.230 


3.216 


1.463 


6399 


4.25 


2.500 


1.366 


0.25 


6400 


3.218 


0.033 


0.632 


0.334 


1.5 


0.00000 


0.0000 


1.233 


3.215 


1.463 


6398 


4.24 


3.500 
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Table 8. Summary of Observational Data for v And 



T e ff (K) [Fe/H] (solar) logg (cgs) vsini (km s 1 ) Reference 

6107 ± 80 0.09 ±0.06 4.01 ±0.10 9.5 ± 0.8 Fuhrmann et al. 1998 

6250 ± 100 0.17 ±0.08 4.30 ±0.10 9.0 ± 0.5 Gonzalez 1997a, 1998 

6125 0.06 3.98 Gratton et al. 1996 

6187 -0.02 4.13 Marsakov et al. 1995 

6205 Blackwell et al. 1994 

6212 0.09 4.17 Edvardsson et al. 1993b 

6198 -0.03 4.22 Balachandran 1990 

6050 0.06 4.0 Boesgaard & Lavery 1986 

9.0 ± 0.4 Gray 1986 

6146 -0.20 4.60 Thevenin & Foy 1983 

6000 -0.23 3.91 Hearnshaw 1974 

6072 -0.14 4.10 Spite & Spite 1973 

6072 -0.11 Powell 1970 
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Table 9. Stellar Models for v And 



M* 




T c ff 


L, 


Z 


X 


Y 


Ago 


M ce 


PSH 




L, 


R, 


T a ff 


log9 


SY 

sz 


Solar 


Solar 


K 


Solar 








Gyr 


Solar 


Solar 


Solar 


Solar 


Solar 


K 


cgs 


Solar 



1.359 


0.17 


6250 


3.333 


0.028 


0.689 


0.283 


2.5 


0.00204 


0.0841 


1.304 


3.336 


1.561 


6252 


4.18 


0.000 


1.347 


0.17 


6250 


3.333 


0.028 


0.676 


0.297 


2.5 


0.00224 


0.0856 


1.297 


3.334 


1.561 


6250 


4.18 


1.500 


1.339 


0.17 


6250 


3.333 


0.028 


0.667 


0.306 


2.6 


0.00238 


0.0866 


1.293 


3.332 


1.561 


6250 


4.18 


2.500 


1.331 


0.17 


6250 


3.333 


0.028 


0.658 


0.315 


2.6 


0.00251 


0.0876 


1.289 


3.330 


1.561 


6249 


4.17 


3.500 


1.353 


0.17 


6150 


3.333 


0.028 


0.689 


0.283 


2.9 


0.00376 


0.0964 


1.310 


3.335 


1.613 


6150 


4.15 


0.000 


1.339 


0.17 


6150 


3.333 


0.028 


0.676 


0.297 


3.0 


0.00406 


0.0979 


1.303 


3.332 


1.613 


6149 


4.15 


1.500 


1.330 


0.17 


6150 


3.333 


0.028 


0.667 


0.306 


3.0 


0.00426 


0.0989 


1.299 


3.330 


1.612 


6148 


4.15 


2.500 


1.320 


0.17 


6150 


3.333 


0.028 


0.658 


0.315 


3.1 


0.00445 


0.0998 


1.294 


3.327 


1.612 


6147 


4.14 


3.500 


1.370 


0.17 


6350 


3.333 


0.028 


0.689 


0.283 


2.0 


0.00000 


0.0000 


1.291 


3.340 


1.513 


6353 


4.22 


0.000 


1.357 


0.17 


6350 


3.333 


0.028 


0.676 


0.297 


2.1 


0.00000 


0.0000 


1.286 


3.337 


1.513 


6352 


4.21 


1.500 


1.349 


0.17 


6350 


3.333 


0.028 


0.667 


0.306 


2.1 


0.00003 


0.0012 


1.283 


3.336 


1.513 


6351 


4.21 


2.500 


1.341 


0.17 


6350 


3.333 


0.028 


0.658 


0.315 


2.1 


0.00031 


0.0144 


1.280 


3.334 


1.513 


6350 


4.21 


3.500 


1.315 


0.09 


6250 


3.333 


0.023 


0.694 


0.283 


2.9 


0.00241 


0.0871 


1.292 


3.330 


1.561 


6248 


4.17 


0.000 


1.310 


0.09 


6250 


3.333 


0.023 


0.687 


0.289 


2.9 


0.00248 


0.0878 


1.290 


3.330 


1.561 


6248 


4.17 


1.500 


1.306 


0.09 


6250 


3.333 


0.023 


0.683 


0.294 


2.9 


0.00253 


0.0882 


1.288 


3.330 


1.561 


6248 


4.17 


2.500 


1.303 


0.09 


6250 


3.333 


0.023 


0.679 


0.298 


2.9 


0.00258 


0.0886 


1.286 


3.329 


1.561 


6248 


4.17 


3.500 


1.403 


0.25 


6250 


3.333 


0.033 


0.684 


0.283 


2.2 


0.00183 


0.0815 


1.313 


3.335 


1.561 


6251 


4.20 


0.000 


1.384 


0.25 


6250 


3.333 


0.033 


0.662 


0.305 


2.2 


0.00207 


0.0837 


1.305 


3.336 


1.561 


6251 


4.19 


1.500 


1.371 


0.25 


6250 


3.333 


0.033 


0.647 


0.320 


2.3 


0.00223 


0.0851 


1.299 


3.337 


1.561 


6251 


4.19 


2.500 


1.358 


0.25 


6250 


3.333 


0.033 


0.632 


0.334 


2.3 


0.00240 


0.0865 


1.293 


3.337 


1.561 


6251 


4.18 


3.500 


1.340 


0.17 


6250 


3.122 


0.028 


0.689 


0.283 


2.5 


0.00218 


0.0821 


1.259 


3.127 


1.511 


6252 


4.21 


0.000 


1.329 


0.17 


6250 


3.122 


0.028 


0.676 


0.297 


2.5 


0.00237 


0.0835 


1.254 


3.125 


1.511 


6251 


4.20 


1.500 


1.321 


0.17 


6250 


3.122 


0.028 


0.667 


0.306 


2.6 


0.00249 


0.0845 


1.250 


3.124 


1.511 


6251 


4.20 


2.500 


1.313 


0.17 


6250 


3.122 


0.028 


0.658 


0.315 


2.6 


0.00261 


0.0854 


1.246 


3.123 


1.511 


6250 


4.20 


3.500 


1.378 


0.17 


6250 


3.558 


0.028 


0.689 


0.283 


2.5 


0.00199 


0.0861 


1.350 


3.563 


1.613 


6252 


4.16 


0.000 


1.367 


0.17 


6250 


3.558 


0.028 


0.676 


0.297 


2.5 


0.00214 


0.0875 


1.344 


3.561 


1.613 


6251 


4.16 


1.500 


1.359 


0.17 


6250 


3.558 


0.028 


0.667 


0.306 


2.5 


0.00224 


0.0885 


1.340 


3.560 


1.613 


6251 


4.16 


2.500 


1.351 


0.17 


6250 


3.558 


0.028 


0.658 


0.315 


2.6 


0.00235 


0.0895 


1.337 


3.559 


1.613 


6250 


4.15 


3.500 
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Tabic 10. Summary of Observational Data for p 1 Cnc 



T e fj (K) [Fe/H] (solar) logg (cgs) usini (km s 1 ) Reference 

5336 ± 90 0.40 ±0.07 4.47 ±0.10 2.5 ±1.0 Fuhrmann et al. 1998 

5250 ± 70 0.45 ± 0.03 4.40 ± 0.15 Gonzalez 1998a 

5150 ± 75 0.29 ± 0.06 4.15 ±0.05 1.4 ±0.5 Gonzalez 1997a 

1 ± 1 Baliunas et al. 1997 

0.414 ± 0.096 Taylor 1991 d 

5100 ± 150 0.20 Arribas & Martinez-Roger 1989 

5350 0.44 Campbell 1978 

5200 0.24 ± 0.09 4.5 Perrin 1977 

5140 0.30 4.4 Oinas 1977 

5196 0.11 4.4 Oinas 1974 

4460 -0.15 Bakos 1971 



d This value is based on a statistical analysis of previous measurements and not any new 
observations. 
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Table 11. Stellar Models for p 1 Cnc 



M. [^1 T off L, Z X Y Ago M ce PSH fl ce L, fl« T off log g |£ 

Solar Solar K Solar Gyr Solar Solar Solar Solar Solar K cgs Solar 



1.029 0.45 5300 0.608 0.053 0.613 0.334 2.4 0.04444 0.0798 0.647 0.612 0.927 5308 4.52 1.500 

0.951 0.45 5300 0.608 0.053 0.579 0.368 8.4 0.04583 0.0825 0.628 0.610 0.927 5305 4.48 2.500 

0.873 0.45 5300 0.608 0.053 0.544 0.403 14.3 0.04722 0.0853 0.609 0.609 0.927 5303 4.45 3.500 

1.048 0.39 5275 0.608 0.046 0.671 0.283 1.7 0.04703 0.0809 0.656 0.606 0.936 5270 4.52 0.000 

0.971 0.39 5275 0.608 0.046 0.630 0.324 7.9 0.04798 0.0835 0.637 0.608 0.936 5275 4.48 1.500 

0.920 0.39 5275 0.608 0.046 0.602 0.351 12.1 0.04862 0.0851 0.624 0.610 0.936 5279 4.46 2.500 

0.868 0.39 5275 0.608 0.046 0.575 0.379 16.3 0.04925 0.0868 0.611 0.611 0.936 5282 4.43 3.500 

0.978 0.45 5325 0.608 0.053 0.579 0.368 5.4 0.04404 0.0803 0.631 0.610 0.919 5329 4.50 2.500 

0.896 0.45 5325 0.608 0.053 0.544 0.403 11.6 0.04466 0.0829 0.613 0.611 0.919 5330 4.47 3.500 

1.006 0.42 5300 0.608 0.049 0.622 0.329 4.1 0.04541 0.0808 0.642 0.609 0.927 5302 4.51 1.500 

0.945 0.42 5300 0.608 0.049 0.591 0.360 8.9 0.04620 0.0828 0.627 0.609 0.927 5303 4.48 2.500 

0.885 0.42 5300 0.608 0.049 0.560 0.390 13.7 0.04699 0.0847 0.612 0.609 0.927 5303 4.45 3.500 

1.045 0.48 5300 0.608 0.057 0.603 0.340 1.2 0.04348 0.0792 0.651 0.616 0.928 5315 4.52 1.500 

0.959 0.48 5300 0.608 0.057 0.565 0.378 7.7 0.04580 0.0823 0.629 0.609 0.927 5303 4.48 2.500 

0.872 0.48 5300 0.608 0.057 0.527 0.416 14.2 0.04811 0.0854 0.607 0.603 0.927 5290 4.45 3.500 

0.991 045 5300 0.569 0.053 0.579 0.368 2Ji 0.04437 0.0775 0.622 0.573 0.897 5308 4~53 2.500 

0.898 0.45 5300 0.569 0.053 0.544 0.403 10.1 0.04585 0.0807 0.601 0.568 0.897 5298 4.49 3.500 

1.020 0.45 5325 0.650 0.053 0.613 0.334 4.7 0.04403 0.0826 0.657 0.653 0.950 5330 4.49 1.500 

0.950 0.45 5325 0.650 0.053 0.579 0.368 9.7 0.04559 0.0850 0.640 0.650 0.950 5325 4.46 2.500 

0.881 0.45 5325 0.650 0.053 0.544 0.403 14.7 0.04714 0.0874 0.622 0.648 0.950 5320 4.43 3.500 
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Table 12. Summary of Observational Data for p CrB 

T e ff (K) [Fe/H] (solar) logg (cgs) -usini (km s _1 ) Reference 

5821 ±80 -0.24 ±0.08 4.12 ±0.10 1.0 ±1.0 Fuhrmann et al. 1998 

5750 ± 75 -0.29 ± 0.06 4.10 ±0.05 ~ 1.5 Gonzalez 1998 

4.27 Kunzl et al. 1997 

5745 -0.22 4.11 Gratton et al. 1996 

5905 -0.24 4.20 Marsakov et al. 1995 

5782 -0.26 4.24 Edvardsson et al. 1993b 

1.5 ± 1.0 Soderblom 1982 

5780 -0.17 3.98 Hearnshaw 1974 

5860 -0.14 Alexander 1967 

5663 -0.20 Wallerstein 1962 
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Table 13. Stellar Models for p CrB 



M, [^] T eff L» Z X Y Ago M ce PSH K ce L. R, T off logg |£ 

Solar Solar K Solar Gyr Solar Solar Solar Solar Solar K cgs Solar 



0.874 -0.23 5750 1.796 0.011 0.706 0.283 14.6 0.03400 0.1256 0.873 1.796 1.354 5750 4.12 0.000 

0.886 -0.23 5750 1.796 0.011 0.718 0.271 14.3 0.03341 0.1250 0.879 1.796 1.354 5750 4.12 1.500 

0.894 -0.23 5750 1.796 0.011 0.725 0.264 14.1 0.03301 0.1247 0.883 1.796 1.354 5750 4.13 2.500 

0.902 -0.23 5750 1.796 0.011 0.733 0.256 14.0 0.03262 0.1243 0.887 1.795 1.354 5750 4.13 3.500 

0.863 -0.23 5675 1.796 0.011 0.706 0.283 15.5 0.04391 0.1334 0.865 1.795 1.390 5675 4.09 0.000 

0.875 -0.23 5675 1.796 0.011 0.718 0.271 15.3 0.04330 0.1330 0.871 1.795 1.390 5674 4.09 1.500 

0.883 -0.23 5675 1.796 0.011 0.725 0.264 15.1 0.04290 0.1326 0.875 1.794 1.390 5674 4.10 2.500 

0.890 -0.23 5675 1.796 0.011 0.733 0.256 15.0 0.04249 0.1323 0.879 1.794 1.390 5673 4.10 3.500 

0.888 -0.23 5825 1.796 0.011 0.706 0.283 13.4 0.02581 0.1177 0.882 1.797 1.319 5825 4.15 0.000 

0.900 -0.23 5825 1.796 0.011 0.718 0.271 13.2 0.02536 0.1171 0.887 1.796 1.319 5825 4.15 1.500 

0.909 -0.23 5825 1.796 0.011 0.725 0.264 13.0 0.02507 0.1167 0.891 1.796 1.319 5825 4.16 2.500 

0.917 -0.23 5825 1.796 0.011 0.733 0.256 12.8 0.02477 0.1163 0.895 1.796 1.319 5825 4.16 3.500 

0.896 -0.17 5750 1.796 0.013 0.704 0.283 13.8 0.03231 0.1242 0.884 1.796 1.354 5750 4.13 0.000 

0.906 -0.17 5750 1.796 0.013 0.713 0.274 13.6 0.03191 0.1237 0.888 1.796 1.354 5750 4.13 1.500 

0.912 -0.17 5750 1.796 0.013 0.720 0.268 13.5 0.03163 0.1234 0.891 1.796 1.354 5750 4.13 2.500 

0.919 -0.17 5750 1.796 0.013 0.726 0.262 13.4 0.03136 0.1231 0.894 1.796 1.354 5750 4.14 3.500 

0.862 -0.23 5750 1.673 0.011 0.706 0.283 15.3 0.03306 0.1208 0.846 1.673 1.307 5750 4.14 0.000 

0.875 -0.23 5750 1.673 0.011 0.718 0.271 15.0 0.03258 0.1203 0.852 1.673 1.307 5750 4.15 1.500 

0.883 -0.23 5750 1.673 0.011 0.725 0.264 14.8 0.03226 0.1200 0.855 1.673 1.307 5750 4.15 2.500 

0.891 -0.23 5750 1.673 0.011 0.733 0.256 14.6 0.03194 0.1197 0.859 1.673 1.307 5750 4.15 3.500 

0.887 -0.23 5750 1.928 0.011 0.706 0.283 13.8 0.03464 0.1303 0.903 1.928 1.403 5750 4.09 0.000 

0.899 -0.23 5750 1.928 0.011 0.718 0.271 13.6 0.03403 0.1297 0.909 1.928 1.403 5750 4.10 1.500 

0.907 -0.23 5750 1.928 0.011 0.725 0.264 13.4 0.03362 0.1293 0.913 1.928 1.403 5750 4.10 2.500 

0.915 -0.23 5750 1.928 0.011 0.733 0.256 13.3 0.03321 0.1289 0.917 1.928 1.403 5750 4.11 3.500 




-10 1 2 3 4 

log 10 (P/days) 



